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ABSTRACT 


This  report  provides  a  brief  summary  of  the  principal  results  obtained  in  a  research 
program  on  the  mechanics  of  progressive  cracking  In  ceramic  matrix  composites  and 
laminates.  The  report  concentrates  on  (1)  progressive  transverse  matrix  cracking  in  cross-ply 
laminates.  (U)  the  effect  of  transverse  matrix  cracks  on  the  axial  response  of  unidirectional 
ceramic  matrix  composites,  (ill)  thermal  conductivities  of  hot  pressed  SiC/BN  composites.  (lv) 
microcracking  in  polycrystalline  ceramics,  and  (v)  the  effect  of  matrix  cracking  and  fiber- 
matrix  interfacial  debonding  on  the  response  of  unidirectional  ceramic  matrix  composites. 
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L  INTRODUCTION 

This  research  project  addressed  several  basic  problems  in  the  damage  mechanics  of 
brittle  composite  materials.  Many  new  results  were  obtained  in  the  course  of  this  work.  The 
principal  results  may  be  summarized  as  follows: 

1.  The  formulation  of  a  complete  model  to  predict  both  first  ply  failure  and  the 
subsequent  progressive  transverse  cracking  in  cross-ply  composite  laminates. 
The  accuracy  of  the  predictions  for  AS-350 1-06  and  T300/934  systems  is 
particularly  gratifying.  Both  theoretically  and  experimentally  it  is  clear  that 
damage  development  in  such  brittle  systems  can  be  retarded  by  keeping  the  ply 
thickness  to  a  minimum. 

2.  The  formulation  of  a  refined  shear  lag  model  for  steady  state  matrix  cracking  in 
unidirectional  ceramic  matrix  composites.  Amongst  other  things,  the  analysis 
shows  the  sensitivity  of  the  model  to  the  mechanics  of  the  interfacial  debond 
regime. 

3.  The  formulation  of  both  the  self  consistent  model  and  the  differential  scheme 
for  the  conductivity  of  three  phase  composites.  In  the  case  of  some  hot  pressed 
SiC/BN  composites  containing  voids,  the  theoretical  predictions  were  compared 
with  experimental  data.  For  the  conductivity  parallel  to  the  hot-pressing 
direction  agreement  was  good.  However  for  the  conductivity  perpendicular  to 
the  hot-pressing  direction  agreement  was  not  good. 

4.  The  investigation  of  the  two-dimensional  hexagonal  array  model  for 
polycrystalline  aggregates.  It  turns  out  that  one  can  obtain  an  analytic  solution 
of  the  problem  for  arbitrary  grain  orientation  distributions.  One  significant 
conclusion  to  emerge  is  that  accurate  determination  of  interfacial  stresses 
demands  consideration  of  at  least  200  contiguous  grains.  In  addition  a  model 
was  formulated  to  describe  progressive  grain  boundary  microcracking.  The 
model  has  the  advantages  of  simplicity,  ease  of  calculation  of  process  zone 
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microcrack  shielding,  allows  for  anisotropy  of  microcracking  and  compares 
favorably  with  models  proposed  by  others. 

5.  The  proper  formulation  and  solution  of  stress  analysis  problems  associated 
with  interfacial  debonding  in  unidirectional  ceramic  matrix  composites.  We 
take  the  opportunity  to  correct  errors  which  are  prevalent  in  the  literature- 
even  In  the  simplest  case  of  an  arc  crack  in  a  homogeneous  material.  The 
solution  to  the  full  problem  is  obtained  from  two  coupled  singular  integral 
equations  and  thus  provides  a  proper  basis  for  all  subsequent  work  on  fiber- 
matrix  Interfacial  debonding. 

These  principal  findings,  together  with  other  related  results,  are  described  In  the  sequel. 


2.  SIGNIFICANT  ACHIEVEMENTS 

2.1  PrtttfrtHuriv*  imiihwm  gjgMflfc ir*  cro— -ply  Gr/cpoxr  and  BJm  cpoanr  laminates. 

In  a  series  of  reports  and  papers  which  were  produced  under  an  earlier  AFOSR  grant  (84- 
0366),  Laws  and  Dvorak  (1,  2,  3)  gave  extensive  results  on  the  Incidence  of  the  first  transverse 
crack  in  cross-ply  laminates  and  the  subsequent  damage  development  under  Increasing  load. 

A  feature  of  the  Laws-Dvorak  analysis  was  that  the  shear  lag  parameter  was  chosen  to  give  the 
correct  value  of  the  stress  for  first  ply  failure.  This  was  entirely  due  to  the  fact  that  the  stress 
Intensity  reduction  factors  required  by  the  Laws-Dvorak  model  [1,2]  were  then  unknown. 

In  the  course  of  the  present  work  this  deficiency  has  been  remedied  by  Laws  and  Wang 
[4],  By  making  use  of  the  analysis  contained  In  Wang's  Ph.D.  dissertation  [5],  Laws  and  Wang  [4] 
have  given  a  completely  deductive  model  for  the  prediction  of  progressive  transverse  cracking 
in  cross-ply  composite  laminates.  It  is  of  interest  to  observe  the  comparison  between  the  Laws- 
Wang  model,  those  of  other  authors  and  the  experimental  data.  Some  results  for  first  ply 
failure  are  shown  in  Figs.  1.  2  and  3.  Also  included  In  Figs.  1,  2  and  3  are  results  obtained  from 
the  theoretical  models  of  Bailey  et  al  [6],  Flaggs  [7],  Nuismer  and  Tan  [8]  and  Hahn,  Han  and 
Croman  [9],  Whilst  the  author  has  reservations  on  the  status  of  the  models  of  these  various 
authors  [6,  7,  8,  9]  it  is  clear  from  Figs.  1,  2  and  3  that  no  model  has  a  decisive  advantage  In  the 
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prediction  of  first  ply  failure--at  least  as  far  as  is  shown  by  the  data  reported  In  these  figures. 

The  situation  Is  entirely  different  when  we  discuss  progressive  transverse  cracking. 
When  we  wish  to  predict  crack  density  as  a  function  of  applied  load  it  is  clear  that  the  models 
given  in  16.  7.  8.  9)  are  not  successful.  Indeed  the  only  acceptable  models  known  to  the  author, 
are  the  Laws -Wang  improvement  of  the  Laws- Dvorak  model  together  with  the  model  of  Wang 
and  Crossman  [11].  It  is  noteworthy  that  predictions  of  both  models  are  almost  identical,  as 
shown  in  Figs.  4  and  5,  wherein  the  models  are  compared  with  data  by  Wang  [1 1]  for  AS-3501- 
06  and  T300/934  systems. 

Finally  we  note  that  a  significant  Implication  of  our  results  Is  that  the  amount  of 
damage  in  laminates  can  be  reduced  bv  keeping  the  plv  thickness  to  a  minimum. 

2£  The  effect  of  transverse  matrix  cracks  on  the  axial  nvmtmac  of  unidirectional  ceramic 

matrix  composites 

The  work  here  has  concentrated  on  an  extension  of  the  Laws-Dvorak  shear  lag  [2] 
model  to  the  steady  state  cracking  regime  discussed  by  Budiansky,  Hutchinson  and  Evans  [  12]. 
We  also  note  that  the  steady  state  cracking  model  has  been  extensively  studies  by  others— most 
notably  by  Dharanl.  Chai  and  Pagano  [13]  and  by  McCartney  [14].  In  this  context  it  is  difficult 
to  compare  the  various  theoretical  models  and  to  compare  with  such  data  as  is  available.  We 
are  therefore  content  to  describe  some  of  the  significant  results  which  emerge  from  the  model 
developed  by  the  author  [15]  and  to  compare  our  results  with  those  of  Budiansky.  Hutchinson 
and  Evans  [12]. 

First  consider  unbonded,  frlctlonally  constrained  fibers.  Then  the  non-dimensional 
critical  cracking  stress  obtained  from  our  model  is  compared  with  the  BHE  model  in  Figs.  6 
and  7.  Aside  from  the  slight  overshoot  in  Fig.  6,  the  two  models  are  comparable.  But 
differences  emerge  when  we  consider  the  critical  matrix  stress  for  Initially  bonded,  debonding 
fibers.  In  this  case  our  model  Is  distinctly  different  from  the  BHE  model.  At  the  expense  of  over 
simplification  it  is  appropriate  to  say  that  our  model  [15]  requires  knowledge  of  the  interfacial 
shear  strength.  Thus  consider  initially  bonded,  debonding  fibers  which  are  then  constrained 
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by  friction.  The  critical  stress  for  various  ratios  of  G<j  Gm  (debonding  toughness/matrix 
toughness)  and  td/xs  (interfacial  shear  strength/friction  stress)  are  shown  in  Figs.  8.  9  and  10. 
In  order  to  examine  the  effect  of  residual  stress  on  the  critical  stress.  Figs.  1 1.  12  and  13  show 
the  different  critical  stresses  when  the  fiber-matrix  interface  in  the  debonded  zone,  remains  in 
contact  with  friction  or  separates.  We  note  that  frictional  contact  produces  a  significant 
increase  in  the  critical  matrix  cracking  stress  as  is  only  to  be  expected. 

2L3  Thermal  conductivities  of  hot  pressed  SIC /BN  re*  tuning 

The  most  significant  results  obtained  in  this  part  of  our  work  has  been  our  success  in 
formulating  the  s.c.m.  and  d.  s.  models  together  with  some  success  in  predicting  the 
experimental  results  of  Ruh,  Bentsen  and  Hasselman  (16)  parallel  to  the  direction  of  hot 
pressing.  However,  as  might  be  expected,  the  theoretical  results  differ  rather  widely  from  the 
experimental  data  for  the  direction  perpendicular  to  the  hot-pressing  direction.  The  results 
shown  in  Figs.  14-17  provide  confirmation  of  the  above  statements. 

Of  course,  a  major  problem  posed  by  any  examination  of  the  data  of  Ruh.  Bentsen  and 
Hasselman  [16]  is  that  they  pertain  to  3-phase  composites.  Indeed,  as  far  as  the  author  is 
aware,  the  work  reported  here  is  the  first  attempt  to  correlate  3-phase  theoretical  models  with 
experiment. 

It  is  abundantly  clear  that  complete  experimental  data  is  essential  for  the  successful 
application  of  any  theoretical  model. 

2.4  MtoncnsMag  In  wlycrartilUnc  ceramist. 

In  the  first  instance  we  mention  the  completion  of  some  earlier  work  (partly  funded  by 
ALCOA)  on  the  effect  of  residual  stress  in  polycrystalline  ceramics.  The  model  proposed  by 
Laws  and  Lee  [17J  is  an  extension  of  the  Evans  [181  two-dimensional  hexagonal  array  model.  It 
is  assumed  that  each  grain  is  elastically  isotropic  but  thermally  anisotropic.  The  orientation 
of  the  various  grains  in  the  array  is  arbitrary. 

It  is  possible  [171  to  give  an  exact  solution  to  this  problem.  Amongst  other  things,  the 


5 


solution  extends  a  result  of  Evans  [18]  to  show  that  the  stress  singularity  at  triple  points  Is 
always  logarithmic.  In  addition,  we  were  able  to  show  that  If  accurate  information  on  the 
residual  stress  at  a  given  interface  is  required,  it  is  essential  to  consider  at  least  2C^ 
surrounding  grains.  This  result  is  in  marked  contrast  to  an  earlier  assertion  by  Evans  and  Fu 
[19.  10.21}. 

A  result  of  significant  practical  interest  is  the  grain  size  for  spontaneous 
microcracking  during  cooldown.  As  described  by  Laws  and  Lee  [17]  the  model  compares 
favorably  with  the  model  of  Hutchinson  and  Tvergaard  [12]  and  the  experimental  data  of  Rice 
and  Pohanka  [23]. 

A  further  area  of  concern  is  the  progressive  microcracking  of  the  polycrystalline 
aggregate  under  continued  mechanical  loading.  At  this  juncture  it  is  perhaps  important  to 
emphasize  that  the  development  of  models  for  progressive  microcracking  of  polycrystalline 
ceramics  play  an  important  role  in  our  work  on  ceramic  matrix  composites. 

The  data  shows  that  such  systems  are  essentially  elastic  in  the  sense  that  residual 
strains  In  simple  tension  tests  are  often  negligible.  Thus  one  is  led  to  formulate  models  for 
microcracking  solids  but  which  have  a  macroscopic  energy  function.  It  is  noteworthy  that 
Hutchinson  [24]  and  Charalambides  and  McMeeklng  [25]  have  emphasized  that  existence  of  a 
microscopic  energy  density  does  not  imply  existence  of  a  macroscopic  energy  density  in  a 
microcracking  solid.  In  fact,  as  is  discussed  by  Laws  [26]  severe  restrictions  must  be  placed  on 
the  microcrack  nucleation  function  in  order  that  the  microcracking  polycrystalline  aggregate 
can  be  regarded  as  macroscopically  hyperelastic.  These  restrictions  give  rise  to  a  nucleation 
function  which  closely  approximates  the  nucleation  function  proposed  by  Charalambides  and 
McMeeklng  [25].  In  effect,  the  different  nucleation  functions  are  demanded  by  the  assumed 
model  for  stiffness  loss  (e.g.  self-consistent  model,  linear  extrapolation,  differential  scheme 
etc).  The  technique  has  considerable  advantages  over  earlier  work  in  that  it  applies  with 
comparative  ease  to  anisotropic  distributions  of  microcracks. 

An  issue  of  much  current  interest  relates  to  the  effect  of  shielding  due  to  microcracking 
in  the  process  zone  of  a  stationary  macroscopic  crack.  This  problem  has  been  discussed  by 
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Hutchinson  (241,  Charalambldes  and  McMeeking  (251  Ortiz  [271  and  others.  An  assessment  and 
comparison  of  the  respective  models  is  discussed  by  the  present  writer  in  [261.  For  the  present 
it  suffices  to  say  that,  within  the  range  of  common  applicability,  the  various  models  are  in 
reasonable  agreement.  The  advantage  which  is  claimed  for  the  writer’s  model  [261,  Is  that  the 
extent  of  microcrack  shielding  is  obtained  with  relatively  little  effort.  In  addition  the  model 
appears  to  have  potential  for  application  to  more  complicated  systems. 

2J5  The  effect  of  matrix  glfiklag  m^r-matrtr  interfadiil  dehondintf  on  the  response  <rf 

linMteSflBMl  matrix  composites. 

It  is  appropriate  to  repeat  in  this  final  report  many  stimulating  conversations  with 
Drs.  Ted  Nicholas  and  Nick  Pagano  and  members  of  their  groups  at  WRDC/MLLN.  In  addition 
it  is  important  to  record  the  excellent  work  by  Zawada  and  Butkus  [281. 

As  stated  in  earlier  reports,  an  extensive  round  table  discussion  was  held  at 
WRDC/MLLN  at  which  it  was  concluded  that  some  exact  solutions  relevant  to  fiber-matrix 
debonding  were  highly  desirable.  Such  solutions  are  extremely  important  for  many  reasons. 

Thus  a  significant  effort  has  been  expended  in  solving,  and  applying,  the  solutions  of  a 
variety  of  fiber-matrix  interfacial  crack  problems.  For  simplicity  I  concentrate  on  the  two- 
dimensional  problem  associated  with  a  single  SiC  fiber  in  an  LAS  matrix  when  there  is  an 
applied  simple  tension  at  infinity  and  there  is  a  single  debond  crack  at  the  interface.  But 
before  I  describe  the  SiC  flber-LAS  matrix,  let  me  briefly  discuss  the  associated  problem 
wherein  fiber  and  matrix  are  of  the  same  material.  The  solution  of  this  problem  is  given  by 
Mushkelishvlli  [291 .  This  solution  was  used  by  Sih,  Paris  and  Erdogan  [30J  to  calculate  the 
stress  Intensity  factors  at  the  crack  tips.  Unfortunately,  as  noted  by  Savin  [31)  the  expression 
for  the  S.I.F.  given  in  [301  is  wrong.  However,  a  significant  feature  of  the  Mushkelishvlli 
solution,  which  appears  to  have  escaped  the  attention  of  all  previous  workers  ,  except  Toya  [321. 
is  that  it  implies  crack  closure  for  certain  orientations  of  the  applied  load,  see  Fig.  18.  This 
unhappy  situation  has  been  remedied  by  Chao  and  Laws  [331  who  gave  an  exact  solution  of  the 
problem  for  partially  closed  cracks.  In  Fig.  19  we  present  a  graphic  to  show  the  extent  of  crack 
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closure  under  changing  orientation  of  applied  load.  Obviously  it  Is  inappropriate  to  give 
comprehensive  details  of  the  techniques  or  results  In  this  report.  It  Is.  however,  instructive  to 
show  the  comparison  between  the  stress  Intensity  factors  given  by  Rooke  and  Cartwright  (34) 
(who  use  the  Savin  [311  result)  and  the  results  obtained  by  Chao  and  Laws  [331.  Despite  the  fact 
that  the  so-called  handbook  solution  is  not  correct  whenever  crack  closure  occurs,  it  is  obvious 
from  Figs.  20-23  that  the  handbook  result  is  remarkable  close  to  the  correct  result.  Further,  in 
Fig.  24  we  show  the  transverse  stiffness  of  a  solid  containing  a  population  of  identical  arc- 

cracks.  with  crack  densities  €  =  0.05  and  €  =  0.01.  it  is  again  clear  that  the  handbook 
solution  is  extremely  accurate  even  when  crack  closure  occurs.  Thus  we  have  the  significant 
(but  anticlimactlc)  result  that  the  Muskhellshvill  solution  (291  Is  entirely  sufficient  for  most 
practical  purposes. 

Turning  now  to  the  genuine  interfacial  crack  problem.  It  is  well  known  that  this 
Involves  an  enormous  increase  in  complexity  [32.  35,  36.  37.  38.  39).  In  addition  recent  work  by 
Rice  [40|  and  Hutchinson  [41]  have  thrown  considerable  light  on  the  problem  of  a  crack  at  the 
Interface  of  two  half  spaces. 

Of  course,  it  Is  clear  that  in  the  case  of  an  interfacial  (arc)  crack  between  an  SIC  fiber 
and  an  LAS  matrix,  we  can  get  crack  closure  by  two  mechanisms:  first  by  the  so-called 
overlapping  surfaces  phenomenon  and  second  because  of  load  orientation.  Nevertheless,  it 
may  be  shown  that  the  problem  can  be  reduced  to  the  solution  of  two  coupled  singular  integral 
equations  which  must  be  solved  numerically.  It  turns  out  that  the  shape  of  the  interfaclal 
crack  with  partial  contact  is  similar  to  that  shown  in  Fig.  19  for  the  arc  crack— but  detailed 
exposition  is  Inappropriate  here.  It  suffices  to  say  that  for  arbitrary  choice  of  fiber  and  matrix 
one  can  calculate  contact  lengths,  stress  intensity  factors,  energy  release  rates,  loss  of 
stiffness,  etc.,  etc.  The  full  details  are  available  In  the  paper  by  Chao  and  Laws[42).  By  way  of 
Illustration  the  actual  contact  lengths  are  shown  In  Figs.  24,  25— note  the  significant  difference 

in  magnitude  of  the  contact  length  8i  (at  the  "open"  end)  compared  with  §2  (at  the  "closed”  end) . 
Also,  typical  stress  intensity  factors  are  shown  In  Figs.  26  and  27.  Perhaps  the  easiest  graphs 
to  Interpret  are  those  showing  loss  of  transverse  stiffness  (Ex)  as  functions  of  material 
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parameters  and  crack  and  load  geometry,  shown  in  Pigs.  28.  29  and  30.  In  this  context,  it  Is  not 
possible  to  assert  that  one  can  use  the  "classical”  overlapping  solution  with  Impunity. 

These  and  other  Issues  are  discussed,  at  length,  by  Chao  and  Laws  (42,  43].  It  is  also 
worth  recording  here  that  significant  progress  has  been  made  by  the  author  In  analyzing  the 
nature  of  the  near-tlp  zones  of  the  interfacial  cracks.  The  analysis  here  shows  that  there  is  a 
definite  "boundary  layer”  effect-similar  to  that  obtained  by  Dundurs  and  Gautesen  [44J. 

Further,  the  analysis  gives  additional  weight  to  Hutchinson’s  [411  proposal  to  set  p  =  0  and 
thus  avoid  the  terms  which  give  rise  to  overlapping  surfaces  at  the  crack  tips.  The  current, 
albeit  tentative,  conclusion  from  our  exact  analysis  Is  that  Hutchinson's  14 1]  proposal  should 
be  entirely  sufficient  for  practical  purposes. 
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Fig.  2  Strength  of  transverse  ply  of  various  T300/934  laminates,  as  predicted,  and 

experimental  data  by  Flaggs  and  Dural  [7].  Thick  ply  result  of  Dvorak  and  Laws  (II 
denoted  by  —  - — . 
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Strength  of  transverse  pi y  of  various  T300/934  laminates,  as  predicted,  and 
ocperttnental  data  by  Wang  [1 11. 
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Fig.  4  Progressive  cracking  ofAS-3501-06  cross- ply  laminates,  as  predicted,  compared  with 
experimental  data  by  Wang  [11]. 
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Fig.  5  Progressive  cracking  of  T300/934  cross-ply  laminates,  as  predicted,  compared  with 
experimental  data  by  Wang  (1 1]. 
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Fig.  7  Matrix  cracking  stress:  unbonded,  frlctlonally  constrained  fibers. 
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Fig.  8  Matrix  cracking  stress:  Initially  bonded,  ftlcUonally  constrained  flbers. 


o 


\o  /  ssaj^s  [«opiJ3 


22 


Fig.  9  Matrix  cracking  stress:  Initially  bonded,  ftlctlonally  constrained  fibers. 
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Fig.  10  Matrix  cracking  stress:  Initially  bonded,  (fictionally  constrained  fibers. 
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Fig.  1 1  Matrix  cracking  stress:  Initially  bonded,  debonding  fibers:  the  effect  of  residual  stress. 
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Fig.  12  Matrix  cracking  stress:  Initially  bonded,  debonding  fibers:  the  effect  of  residual  stress. 
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Fig.  13  Matrix  cracking  stress:  Initially  bonded,  debonding  fibers:  the  effect  of  residual  stress. 
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Fig.  14  Thermal  conductivity  of  SIC /BN  composites  (containing  voids  of  aspect  ratio  5) 

parallel  to  the  hot-pressing  direction  by  the  self  consistent  method. 
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Fig.  15  Thermal  conductivity  of  SiC/BN  composites  (containing  voids  of  aspect  ratio  5) 
parallel  to  the  hot-pressing  direction  by  the  differential  scheme. 


Volume  Concentration  of  Boron  Nitride 


17  Thermal  conductivity  of  SiC/BN  composites  (containing  voids  of  aspect  ratio  5) 
parallel  to  the  hot-pressing  direction  by  the  differential  scheme. 
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Fig.  19  Shape  of  partially  closed  crack  for  various  load  orientation  with  <j>  ■  60°  (not  to  scale). 
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Fig.  20  Comparison  of  Hi  (A)  for  exact  solution  and  handbook  solution  1341  when  $  » 


Fig.  21  Comparison  of  Kg  (B)  for  exact  solution  and  handbook  solution  [34 1  when  $ 


Fig.  23  Comparison  of  Ka  (B)  for  exact  solution  and  handbook  solution  I34|  when  $  =  45°. 
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Fig.  24  Reduction  in  Young's  modulus  for  given  arc  crack  density  (e).  when  $  =*  60°  comparison 
of  exact  solution  and  handbook  solution  (34|. 
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,  and  load  orientation,  ($,  when  the  first  Dundurs  parameter  lsO*a  0.5  and  $ 
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when  the  first  Dundurs  parameter  a*  ■  0.5  and  $ 
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Fig.  28  Mode  II  S.I.F.  at  B  for  various  values  of  the  Dundurs  parameter,  p*.  and  load 
orientation,  ay  when  the  first  Dundurs  parameter  ol*  =*  0.5  and  $  a  60°. 
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0.75,  $  =  60°  and  the  volume  fraction  of  fibers  ■  0.05. 
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1.  INTRODUCTION 

The  development  OF  a  satisfactory  theory  for  cross-ply  laminates  which 
have  been  damaged  by  transverse  matrix  cracking  under  monotonic  loading 
has  attracted  a  substantial  number  of  investigators.  The  formulation  of  a  shear  lag 
model  appears  to  have  been  first  proposed  in  a  series  of  papers  by  Bailey  and  his 
co-workers  [1,2, 3,4,5$.  This  work,  in  turn,  relies  on  some  studies  of  unidirec¬ 
tional  composites  by  Aveston  and  Kelly  [6].  Subsequent  contributions  to  the 
theory  have  geen  given  by  Wang  [7],  Highsmith  and  Reifcnider  [8],  Flaggs  and 
Kural  [9],  Nuismer  and  Tan  [10],  Manders.  Chou,  Jones  and  Rock  [11], 
Fukunaga,  Chou,  Peters  and  Schulte  [12],  Flaggs  [13],  Ohira  [14]  and  Ogin. 
Smith  and  Beaumont  [13,16].  Doubtless  a  diligent  search  of  the  literature  would 
disclose  other  related  work. 

In  an  important  series  of  papers  Wang  and  Crossman  [17,18]  and  Wang  and  his 
co-workers  [19,20,21]  have  discussed  transverse  cracking  from  a  different  point  of 
view.  And,  importantly  for  the  present  investigation,  the  work  of  these  authors 
contains  some  comprehensive  experimental  data  on  the  progressive  transverse 
cracking  of  composite  laminates. 

Additional  work  on  the  loss  of  stiffness  of  cracked  composite  laminates  which 
is  based  on  three-dimensional  stress  analysis  has  been  given  by  Laws  and  Dvorak 
[22,23,24,25]  and  Hashin  [26]. 

It  is  not  our  aim  here  to  give  a  critical  survey  of  the  existing  literature.  Nor  do 
we  attempt  to  point  out  the  various  similarities  and  differences  between  the  pub¬ 
lished  work  and  the  work  described  below.  Nevertheless,  there  are  some  major 
differences  which  should  be  emphasized  here. 
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In  the  first  place,  the  existence  of  residual  stresses  is  fully  accounted  for 
throughout  the  whole  of  this  analysis.  Actually,  we  show  that  such  stresses 
demand  that  there  is  a  permanent  strain  when  the  applied  load  is  large  enough  to 
cause  transverse  cracking.  Physically  this  result  is  obvious.  However,  it  turns  out, 
both  theoretically  and  experimentally,  that  this  permanent  strain  is  negligible. 
But  perhaps  the  main  difference  between  the  model  proposed  here  and  existing 
models,  lies  in  our  treatment  of  the  statistics  of  progressive  cracking.  In  par¬ 
ticular  we  note  that  Bailey  et  al.  [1-5]  assume  that  cracks  always  occur  midway 
between  existing  cracks,  whereas  Manders  et  al.  [11]  and  Fukunaga  et  al.  [12]  use 
a  Weibull  strength  distribution  of  the  transverse  ply  in  their  discussion  of  progres¬ 
sive  cracking.  On  the  other  hand  Wang  and  Crossman  [17]  introduce  distributions 
of  effective  flaw  sizes  and  locations  in  the  transverse  ply.  By  way  of  contrast,  the 
approach  adopted  here  supposes  that  a  transverse  crack  will  propagate  when  it  is 
energetically  favorable  and  that  the  location  of  this  transverse  crack  is  associated 
with  a  probability  density  function.  Clearly,  a  crucial  issue  is  the  specification  of 
this  probability  density  function.  Based  on  simple  fracture  mechanics  we  suggest 
a  precise  choice  for  the  required  probability  density,  namely  that  it  is  propor¬ 
tional  to  the  stress  in  the  transverse  ply.  Of  course,  normalization  gives  the  factor 
of  proportionality.  This  choice  is  explored  and  shown  to  give  good  agreement 
with  experiment. 

Given  this  choice  of  probability,  the  only  parameter  used  here  which  is  not  de¬ 
termined  by  standard  tests  is  the  so-called  shear  lag  parameter.  However,  we 
show  that  one  can  infer  this  parameter  horn  standard  data  and  knowledge  of  the 
first  ply  failure  stress.  We  therefore,  propose  an  explicit  formula  for  the  deter¬ 
mination  of  the  shear  lag  parameter  £.  We  remark  that  the  values  thus  obtained 
do  not  agree  with  those  obtained  from  a  formula,  due  to  Garrett  and  Bailey  [l]. 

Thus,  the  present  analysis  provides  a  well-defined  model  for  transverse  crack¬ 
ing  in  cross-ply  composite  laminates  based  on  statistical  fracture  mechanics.  This 
model  is  well-defined  in  the  sense  that  no  adjustable  parameters  are  available  to 
fit  a  particular  set  of  experiments.  Amongst  other  things,  the  model  delivers 
explicit  formulae  for  the  loss  of  stiffness  as  a  function  of  crack  density,  and  for 
crack  density  as  a  function  of  applied  load. 

2.  BASIC  EQUATIONS 

For  simplicity,  we  are  here  concerned  with  a  strictly  one-dimensional  theory 
of  symmetric  cross-ply  composite  laminates.  In  addition,  we  only  consider 
monotonic  simple  tensile  loading,  see  Figure  1.  Generalization  to  angle-ply  lami¬ 
nates  and  biaxial  loading  will  be  reported  in  a  further  paper. 

It  is  well  known  that  the  strength  of  composite  laminates  depends  upon  the 
residual  stresses  due  to  cool-down.  Thus  let  the  initial  stresses  in  the  laminate  be 
of  and  of  in  the  transverse  and  longitudinal  plies  respectively.  Here,  and  subse¬ 
quently,  the  subscripts  t  and  t  are  used  to  denote  transverse  and  longitudinal 
respectively.  Obviously 


bof  +  dof  =  0 


(1) 
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*b 


where  b  is  the  thickness  of  the  outer  0°  plies  and  d  is  the  half-thickness  of  the 
inner  90°  plies,  see  Figure  1. 

As  usual,  it  is  convenient  to  measure  displacements  from  the  state  of  initial 
stress  with  no  mechanical  loading.  In  shear  lag  theory,  it  is  assumed  that  the  dis¬ 
placement  of  a  particular  layer  is  uniform  over  the  thickness  of  that  layer:  let  u(jr) 
be  the  displacement  of  the  0°  layer  and  v(x)  be  the  displacement  of  the  90°  layer. 
The  associated  strains  of  the  longitudinal  and  transverse  plies  are: 


€, 


du 

dx 


(2) 


Young's  modulus  for  the  uncracked  laminate  is,  in  this  approximation. 


bE,  +  dE, 
**  b  +  d 


(3) 


Let  <r„  a,  be  the  total  stress  in  the  respective  plies,  then  overall  equilibrium  of 
the  laminate  demands  that 


bo,  +  do,  =  (b  +  d)o. 


W 


figure  2.  Stresses  on  individual  layers  in  the  laminate. 
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where  <x.  is  the  applied  stress.  Also  if  r  is  the  shear  stress  in  the  shear  layer. 
Figure  2,  then  equilibrium  of  the  0s  and  90°  plies  separately  implies  that 


do,  da, 

r  =  -b—  =  d  — 
dx  dx 


(5) 


In  this  paper  we  examine  cross  ply  laminates  which  are  subject  only  to 
mechanical  loading  a.  -  as  well  as  residual  stress.  There  is  no  difficulty  in  extend¬ 
ing  the  analysis  to  encompass  thermal  loading  and  this  will  be  reported  else¬ 
where.  The  stress-strain  relations  for  the  respective  plies  are,  therefore. 


a,  =  <rf  +  r,. 


r,  =  £,e. 


a,  =  a f  +  r„ 


r,  =  E,(, 


(6) 


In  other  words,  r,  and  r,  are  the  stresses  due  to  mechanical  loading.  Finally,  the 
essential  ingredient  of  shear  lag  theory  is  that  the  shear  stress  r  is  assumed  to  be 
proportional  to  the  relative  displacement  (v  —  u): 

t  =  K(v  -  u)  (7) 


where  AT  is  a  constant. 

Analysis  of  the  fracture  mechanics  of  transverse  plies  requires  knowledge  of 
the  strain  energy  of  the  laminate -or  more  precisely  the  total  energy  of  the  lami¬ 
nate.  As  a  prelude  to  this  calculation,  we  here  evaluate  the  strain  energy,  per  unit 
width  of  the  laminate,  between  two  sections  PP  and  QQ,  see  Figure  3.  Now  the 
increase  in  strain  energy,  W,  due  to  the  application  of  mechanical  loads  r,  and  r, 
[see  Equation  (4)]  is  equal  to  the  work  done: 


W  =  b[(r,  +  2  o*)u\9  +  d[(r,  +  2  a?)v]? 
=  b[(a,  +  o,*)u)3  +  d[(o,  +  Ov]P 


(8) 


Figure  3.  Section  of  uncrackod  laminate. 
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and  we  have  used  the  notation 

[«]?  =  u(Q )  -  u(P) 

In  order  to  obtain  Equation  (8)  it  is  essential  to  remember  that  the  applied  loads 
on  the  respective  plies  increase  from  zero  to  r,  and  r,,  whereas,  the  residual 
stresses  are  essentially  dead  loads— hence  the  factor  2.  Clearly  we  can  rewrite  (8) 
in  the  form 


d 


((<r,  +  of)u]dx  +  d 


1> 


[(a,  +  trf)v)dx 


With  the  help  of  (2),  (5),  (6),  and  (7)  we  may  show  that 


It  is  now  a  trivial  exercise  to  derive  a  complete  set  of  differential  equations  for 
shear  lag  theory.  In  order  to  discuss  transverse  cracking,  the  most  useful  equation 
governs  the  stress  in  the  transverse  ply.  One  form  of  this  equation  is 


d'o,  11 

dx2  d2 


where  the  non-dimensional  shear  lag  parameter  £  is  given  by 


KdjbE,  -)-  dE,) 
bE,E, 


(10) 


(11) 


Let  us  now  consider  the  straightforward  problem  associated  with  two  trans¬ 
verse  cracks  distant  2 h  apart,  see  Figure  4.  We  are  particularly  interested  in  the 
elastic  field  in  that  part  of  the  laminate  which  is  between  the  cracks.  Thus,  the  ap¬ 
propriate  boundary  conditions  for  Equation  (10)  are 


a,  =  0  when  x  -  ±  h 


The  required  solution  is 


(12) 


(13) 


Progressive  Transverse  Cracking  In  Composite  Laminates 


905 


- - 2d  10 - - 

Figure  4.  Two  adjacent  transverse  cracks  in  the  90 °  plies. 


The  stress  in  the  longitudinal  plies  is 


E ,  „ 

o,  =  T  *.U  + 


. 

dE.  c°shT 

bE,  (h 
cosh  —r 
a 


+  a<  1  ~ 


*  & 
cosh  —r 
d 

.  9i 
cosh  — 
a 


(14) 


From  (2),  (6),  (13)  and  (14)  we  find  that  the  associated  displacements  are 


“  =  Eo  X  +  &E, 


■  h  & 

f°f  *  6  '•> - £ 

cosh  — 
a 


+  cx 


(15) 


a.  d 

V  =  E0X  ~  «£, 


(fff  +  ~rr  a.) 


sinh 


£x 


Ee  .  # 
cosh  — 
a 


+  c. 


(16) 


where  c,  and  c»  are  constants. 


3.  LOSS  OF  STIFFNESS 

With  the  help  of  the  results  of  the  preceding  section  we  can  now  calculate  the 
loss  of  stiffness  of  the  cross-ply  laminate  due  to  transverse  cracking.  Thus,  let  the 
average  distance  between  transverse  cracks  be  2 h.  Figure  4.  Then  the  average 
strain,  e. ,  of  the  uncracked  ligament  AB,  as  measured  at  the  surface  of  the  lami¬ 
nate,  is  given  by 


u(B)  -  u(A) 

2 h  ~  Eq 


d'E,  & 
fr bE,  tanh  d 


d'o* 

&hE, 


(17) 
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where  we  have  used  (15).  We  note  that  we  would  still  recover  (17)  if  we  were  to 
consider  any  two  sections  distant  2 h  apart.  A  trivial  rearrangement  of  (17) 
together  with  use  of  the  crack  density  parameter 


0  =  dlh 


now  yields 


da?  -  )  (,  ^  , 

-bE,  ^ll1  +I*£Ttanh 


0  da? 
* 


(18) 


which  is  the  effective  stress-strain  relation  for  the  cracked  laminate.  We  observe 
that  Equation  (18)  shows  that,  under  cracking,  the  laminate  acquires  a  permanent 
strain,  e,,  due  to  initial  stress: 


e.  = 


0  daf  £ 

1 6r**i 


(19) 


It  would  be  misleading  to  give  the  impression  that  the  permanent  strain  (19) 
were  important  either  experimentally  or  theoretically.  In  fact  for  all  laminates 
considered  here  the  shear  lag  parameter  lies  between  0.5  and  2.5— as  is  discussed 
later.  Accordingly  the  value  of  e,  is  at  most  5  %  of  the  value  of  e.  in  the  worst  pos¬ 
sible  situation,  i.e.  0  large.  Hence  we  can  neglect  e,  in  applications  and  arrive  at 
the  elegant  formula  for  the  Young’s  modulus  E(ff)  of  a  laminate  containing  trans¬ 
verse  cracks  of  density  0: 

Ws  note  that  as  the  crack  density  0  —  0,  E  —  Eo  as  it  should.  In  addition  as 
0  —  oo  we  see  that 


-&(i  + 


bE, 

b  +  d 


which  is  the  result  given  by  ply  discount  theory. 

An  assessment  of  the  validity  of  Equation  (20)  will  be  given  in  Section  5  when 
we  have  evaluated  the  shear  lag  parameter  £. 


4.  PROGRESSIVE  CRACKING 

We  now  turn  to  the  problem  of  determining  the  transverse  crack  density  0  as 
a  function  of  the  increasing  applied  load  <r. .  Consider  the  uncracked  ligament  AB 
as  in  Figure  4.  When  the  applied  load  a.  reaches  a  critical  value,  this  ligament 
will  itself  crack  at  some  location  C  as  in  Figure  5.  There  is  no  reason  to  suppose 
that  C  lies  at  the  mid-point  of  AB.  Assuming  that  the  additional  cracking  is  pro- 
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AC  B 


2h 


Figure  5.  The  ligament  AB  with  an  additional  crack  at  C. 


duced  under  fixed  loads,  we  can  use  the  analysis  of  Section  2  to  calculate  the 
required  energy  release  rate. 

Let  us  temporarily  regard  the  configuration  of  Figure  4  as  state  1  in  which  the 
displacements  are  ut  and  v,,  whilst  the  configuration  of  Figure  5  is  state  2  with 
displacements  u2  and  v2.  From  (8)  the  strain  energy,  per  unit  width,  of  the 
ligament  AB  in  state  1  is: 


Wi  *=b{(o,  +  <r?)ut]J!  +  d\(o,  +  uf)v,]J 


(21) 

=  \(b  +  d)a.  +  bon  MS  +  do ?  [v,]- 

since  a,(A)  =  o,(B)  —  (X  Also  the  combined  strain  energy,  per  unit  width,  of 
the  two  ligaments  AC  and  CB  in  state  2  is,  from  (21), 

W*  =  ((6  +  d)o,  +  bo?]  IMS  +  M?)  +  dfff{[v2]S  +  [v,]?)  (22) 

But  Ui  is  continuous  at  C,  hence 


MS  +  [«,]?  -  MS  (23) 

Under  the  assumption  of  cracking  at  fixed  loads,  the  work  done,  per  unit  width, 
on  the  segment  AB  when  the  extra  crack  occurs  at  C  is 

2b[o,(uj  -  «t)]S  =  2(b  +  d)o.[ut  -  k,]S  (24) 

Thus  the  total  energy,  <5,  per  unit  width,  released  by  introduction  of  the  crack  at 
C  is  the  work  done  by  the  applied  loads  minus  the  increase  of  strain  energy.  From 
(21)  to  (24)  we  have: 


<2=  2(b  +  d)o.[ut  -  «,]2  -  (Wj  -  W,) 

=  [( b  +  d)o.  +  do?)[ui  -  «,]S  +  d<7?[v,  -  v2]S  +  <fa?[v,  -  v2]f 
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A  series  of  routine,  but  tedious,  calculations  shows  that 


2 d\b  ±  d)Eo  ,  ,  £  ' 

&e,e,  (<t'  +  E„ ffm) ; 


&,  ^  &h  & 

tanh  ~2d  +XZnhU  ~  anh  T 


But  the  energy  release  rate  G ,  per  unit  length,  for  crack  propagation  across  the 
ply  is  determined  from  G  -  <SHd,  see  Dvorak  and  Laws  [24],  Thus 


G  m 


d(b  +  d)E0 
&E,E, 


(fff  + 


Z'->‘ 


tanh 


2d 


+  tanh 


Jfh 

2d 


—  tanh 

a 


(25) 


It  is  now  easy  to  obtain  a  first  ply  Mure  criterion  from  (25)  by  taking  A,  A,, 
A,,  to  be  large  and  setting  G  -  Gc.  The  result  is 


&E,EoGc  ]•'*  _  Eo 
d(b  +  d)E,  E. 


(26) 


When  we  are  given  standard  data,  namely  b,  d,  E„  E„  the  coefficients  of  ex¬ 
pansion  and  the  temperature  drop  during  cool-down,  we  easily  calculate  a"  and 
Eo.  Thus  (26)  provides  a  relationship  between  a1!1,  G.  and  £.  However,  a™  and 
Gc  are  readily  measured,  and  we  therefore  regard  Equation  (26)  as  the  rule  which 
determines  the  shear  lag  parameter  £. 

Once  £  has  been  so  determined.  Equation  (20)  provides  an  explicit  formula  for 
the  loss  of  stiffness  of  the  laminate. 

As  for  subsequent  cracking,  a  major  complication  arises  because  the  location 
of  the  ■‘next”  crack  cannot  be  obtained  by  deterministic  methods.  Indeed,  this  is 
precisely  the  situation  indicated  by  experiment.  Thus  we  must  proceed  on  a 
statistical  basis. 

In  order  to  analyze  progressive  cracking,  suppose  that  the  laminate  contains 
transverse  cracks  with  average  separation  2A  with  associated  crack  density 
0  =  d/h  as  in  Figure  4.  Then  the  next  crack  which  appears  in  the  ligament  AB 
will  be  at  some  location  C  as  in  Figure  5.  Let  <x.(A,)  be  the  applied  stress  which 
is  needed  to  produce  the  crack  at  C.  We  can  determine  <r.(A,)  directly  from  (25) 
by  setting  G  —  G,.  When  use  is  also  made  of  (26),  it  follows  that 


<r.(A,)  =  ( a +  £-  a?) 


tanh 


2d 


+  tznh 


lb. 

2d 


~  tanh  “j 


E, 


(27) 


Since  the  location  of  C  is,  in  an  appropriate  sense,  random,  let  p  be  the  prob¬ 
ability  density  function  for  the  site  of  the  next  crack.  Thus  in  a  laminate  which 
already  contains  cracks  of  density  0,  the  expected  value  of  the  applied  stress  to 
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=  j  p(y)a.(y)dy  (28) 

where  o.(y)  is  given  by  (27).  The  choice  of  probability  density  function  is 
crucial.  Three  candidates  immediately  present  themselves. 

Case  1 

An  extreme  situation  would  arise  if  the  next  crack  were  guaranteed  to  occur  at 
the  mid-point  of  the  unbroken  ligament.  This  would  imply  the  choice 

P(y)  -  H  v  -  A)  (29) 

where  5(  y)  is  the  Dirac  delta  function. 

Case  2 

Another  extreme  would  be  to  assume  that  all  locations  in  the  ligament  were 
equally  likely,  giving 

p(y)  =  ^  (30) 

Case  3 

A  more  appealing  hypothesis,  based  on  simple  fracture  mechanics,  would  be 
to  assume  that p{y)  is  proportional  to  the  stress  in  the  transverse  ply: 

ply)  *  <r.(y) 

Thus  from  Equation  (13)  with  x  =  y  —  h 


For  Case  1,  we  can  obtain  the  explicit  solution 

£•  £  t  £■ 

E[a.(0)]  =  {a1:1  +  Y  a' )  2  tanh  jq  ~  ^  ~jf  ~  ~E  (32) 

But  for  Cases  2  and  3,  the  integral  must  be  evaluated  numerically. 
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5.  COMPARISON  WITH  EXPERIMENT 

In  the  first  place  we  compare  the  predictions  of  the  foregoing  theory  with  some 
experimental  data  of  Highs  mi  th  and  Reifsnider  [8].  These  authors  reported  data, 
and  some  theoretical  results,  for  several  E-glass -epoxy  systems.  Of  particular 
concern  here  are  the  data  for  (0,  90j),  specimens. 

The  data  for  the  E-glass-epoxy  systems  studied  in  [8]  are  as  follows 

E,  =  41.7  GPa  E,  =  13.0  GPa 

<rf  =  8.4  MPa  aV  -  55  MPa 

ply  thickness  0.20  mm. 

Unfortunately  Highsmith  and  Reifsnider  [8]  do  not  give  a  value  for  G..  However, 
the  data  given  above  has  been  used  by  Laws  and  Dvorak  [25]  and  by  Hashin  [26] 
to  validate  the  loss  of  stiffness  given  by  the  self-consistent  model  and  a  lower 
bound  prediction  respectively.  The  Hashin  [26]  lower  bound  and  the  experimen¬ 
tal  data  are  displayed  in  Figure  6.  We  remark  that  the  self-consistent  prediction 
is  not  indicated  in  Figure  6  because  it  is  only  marginally  greater  than  the  Hashin 
[26]  bound  and  thus  the  two  curves  are  virtually  indistinguishable. 

But  in  order  to  obtain  the  shear  lag  prediction  for  loss  of  stiffness,  lack  of  a 
definite  value  for  G.  poses  a  problem.  In  this  connection  we  note  that  in  an  earlier 


Figure  8.  Experimental  and  theoretical  values  for  stiffness  loss  of  (0,90JS  E-glass  epoxy ; 

laminate:  (1)  Highsmith-Reitsnidar  prediction . .  (2)  shear  lag - .  (3) 

lower  bound - .  Experimental  data  from  Reference  [BJ. 
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paper  Dvorak  and  Laws  [24]  studied  some  E-glass-epoxy  systems  and  used  the 
value  Gc  =  250  Jm~\  If  we  were  to  use  the  same  value  here,  then  Equation  (26) 
would  imply  the  value  £  =  0.7.  The  loss  of  stiffness  predicted  by  Equation  (20) 
is  indicated  in  Figure  10.  We  note  that  the  resulting  curve  gives  reasonable  agree¬ 
ment  with  experiment  but  violates  the  Hashin  lower  bound  for  small  crack  den¬ 
sities.  We  might  be  prepared  to  accept  this  state  of  affairs  as  being  an  unfortunate 
but  unavoidable  consequence  of  the  approximate  shear  lag  theory.  However,  a 
more  reasonable  interpretation  is  that  the  ad  hoc  choice  of  Gc  —  250  Jm'2  is  not 
correct.  In  order  to  develop  this  line  of  reasoning  we  note  that  the  maximum 
value  of  Gt  which  implies  that  the  predicted  loss  of  stiffness  is  entirely  consistent 
with  the  Hashin  bound  is  Gc  =  193  Jnr2.  This  in  turn  implies  that  £  =  0.9  and 
thus  a  predicted  loss  of  stiffness  which  is  shown  in  Figure  6. 

For  completeness  Figure  6  also  shows  the  comparison  between  the  theoretical 
prediction  of  Highsmith  and  Reifsnider  [8]  with  their  experimental  data. 
Although  these  authors  also  use  a  form  of  shear  lag  theory,  it  is  not  easy  to  com¬ 
pare  their  model  with  the  theory  given  here. 

We  note  that,  for  the  graphite-epoxy  systems  discussed  later,  data  is  complete 
and  there  is  no  ambiguity  in  the  value  of  G.  and  hence  of  £. 

Despite  the  incompleteness  of  the  data  for  the  E-glass  systems,  comparison  of 
theory  with  experiment  for  the  loss  of  stiffness  of  (0,  90]),  laminates  ‘■annot  be 
regarded  as  decisive.  In  the  first  place  the  loss  of  stiffness  of  typical  graphite- 
epoxy  systems  due  to  transverse  micro-cracking  is  minimal.  In  the  second,  the 
major  advantage  of  the  shear  lag  model  lies  in  its  ability  to  provide  an  explicit 
prediction  of  crack  density  as  a  function  of  applied  load.  In  feet,  if  we  agree  to 
determine  £  from  Equation  (26),  shear  lag  theory  does  not  contain  any  adjustable 
parameters.  Rather,  we  need  to  identify  the  conect  probability  density  function. 
Once  the  choice  of  p(y)  has  been  made,  the  statistical  theory  presented  here  is 
well-defined. 

Let  us  now  turn  to  the  Highsmith-Reifsnider  [8]  data  for  crack  density  as  a 
function  of  applied  load.  In  view  of  our  earlier  discussion  we  here  take  G.  =  193 
Jm-1  which  implies  that  £  =  0.9.  Theoretical  results  for  crack  density  as  a  func¬ 
tion  of  applied  load  can  be  obtained  from  Equation  (28)  for  the  three  choices  of 
probability  density  function  (29),  (30)  and  (31).  The  various  curves  and  the  data 
are  shown  in  Figure  7.  We  note  in  passing  that  the  numerical  evaluation  of  some 
of  the  integrals  requires  considerable  attention  to  detail.  Clearly  the  most  promis¬ 
ing  choice  of  probability  function  is  Case  3.  Indeed  one  can  make  a  strong  a 
priori  argument  for  Case  3  based  on  elementary  fracture  mechanics.  We  there¬ 
fore  propose  the  use  of  the  probability  density  function  (31)  as  appropriate  for  the 
prediction  of  progressive  cracking. 

It  is  of  interest  to  show  the  sensitivity  of  the  predicted  crack  density  to  the  value 
of  G. .  This  is  indicated  in  Figure  8  wherein  it  is  clear  that  the  theory  makes  the 
(necessary)  prediction  that  the  tougher  the  material  the  less  the  crack  density  for 
given  load. 

We  now  turn  to  some  different  work  by  Wang  and  Crossman  [17-21].’  The  work 
of  these  authors  is  exclusively  concerned  with  graphite-epoxy  systems  and  is 
centered  on  the  prediction  of  crack  density  as  a  function  of  the  applied  (mono- 


"  ■  -  _  _ _  _  ■  ■  __  »  — / 

**n*i~Om  from  Referee 

[8J.  Predictions  obtained  from  probabtity  distribution  3  for  indicated  values  of  Gc. 
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tonic)  load.  Those  parts  of  the  work  of  these  latter  authors  which  are  devoted  to 
fatigue  are  not  relevant  to  the  present  discussion. 

For  our  present  purpose  the  most  convenient  source  for  data  is  the  survey 
article  fay  Wang  [20].  We  first  consider  the  AS-3501-06  graphite  epoxy  systems. 
Results  for  (0,,  90),,  (0,.  90j),  and  (0},  90,),  laminates  are  given  in  Figures  11, 
12  and  O  of  Wang's  [20]  article.  When  we  make  use  of  the  data  quoted  by 
Wang  [20],  see  also  [21],  we  find  the  values  of  the  shear  lag  parameters  are  as 
follows 

(0».  90),  $  =  0.93 

(0*.  90,),  $  =  1.38 

(0„  90,).  $  =  2.24 

Comparison  of  the  respective  shear  lag  predictions  with  the  experimental  results 
are  shown  in  Figure  9.  The  reader’s  attention  is  drawn  to  the  foot  that  we  have 
omitted  Wuig’s  [20]  numerical  results  from  Figure  9  because  it  is  impossible  to 
do  justice  by  replicating  the  published  graphs.  The  important  observation  is  that 
both  shear  lag  and  the  Wang-Cross  man  theory  give  very  good  predictions.  This 


Applied  load  (MPa) 

Figure  9.  Theory  versus  experiment  for  progressive  cracking  of  AS-3501-06  laminates.  Data 
from  Wang  [20]. 
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Figure  10.  Theory  varans  axporirrmrt  for  progressive  cracking  ot  T300/934  laminatas.  Data 
from  Wang  ( 20/ . 


is  quite  remarkable  since  the  respective  theories  are  based  on  entirely  different 
premisses. 

Finally,  we  turn  to  the  data  reported  by  Wang  [20]  for  some  T300/934  lami¬ 
nates.  Again  we  use  the  data  given  by  Wang  [20,21]  to  obtain  the  following  values 
for  [  for  the  indicated  lay-ups: 


(0,  90,.  0) 

£  =  1.08 

s 

P 

3 

3 

£  =  1.70 

(0,  904,  0) 

£  =  1.79 

The  theoretical  predictions  are  compared  with  the  experimental  data  in  Figure  10. 
Again  it  is  encouraging  to  report  excellent  agreement. 
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Abstract 

A  hexagonal  grain  array  model  is  used  to  study  grain  boundary  microcracking  of  a  polycrystalline 
aggregate  due  to  residual  stress.  Each  grain  is  assumed  to  be  elastically  isotropic  but  thermally  anisotropic. 
The  axes  of  thermal  anisotropy  for  each  grain  are  arbitrary.  An  explicit  analytic  solution  is  obtained  for 
the  entire  residual  stress  field.  This  solution  is  used  to  give  a  detailed  description  of  the  grain  boundary 
stress  fields.  Further,  explicit  algebraic  formulae  are  given  for  stress  intensity  factors  associated  with  grain 
boundary  microcracks.  The  results  are  used  to  predict  the  critical  grain  size  for  the  occurrence  of  spon¬ 
taneous  microcracking.  Agreement  between  theory  and  experiment  is  good. 


1.  Introduction 

This  paper  is  devoted  to  the  study  of  grain  boundary  microcracking  in  single  phase 
polycrystalline  brittle  solids  due  to  cooldown.  It  has  long  been  held  (Clarke.  1964; 
Cleveland  and  Bradt,  1978;  Davtdge,  1981 ;  Evans,  1978;  Fu  and  Evans.  1982, 
1985;  Kuszyk  and  Bradt,  1973;  Ohya  et  al.,  1987)  that  residual  stresses  in 
polycrystalline  aggregates  are  due  to  thermal  expansion  anisotropy — it  being 
invariably  assumed  that  each  grain  is  elastically  isotropic.  There  is  complete  agreement 
that  there  is  a  critical  minimum  grain  size  for  microcracking  to  occur.  Further,  there 
is  agreement  that  for  given  cooldown  temperature  change  (AT),  Young's  modulus  £, 
difference  between  the  larger  and  the  mean  coefficient  of  thermal  expansion  (Aa)  and 
grain  boundary  toughness  (Gfb),  the  formula  for  the  critical  facet  length  (4)  is  of  the 
form 


4  =  2 


£(AaAT):’ 


(1) 


where  Q  is  a  constant.  It  is  not  our  intention  to  give  a  critical  review  of  the  various 
arguments  which  have  been  used  to  arrive  at  proposed  values  for  Q.  All  that  need  be 
said  Is  that  some  authors  use  approximate  stress  analysis  whereas  others  choose  Q  to 
fit  certain  experimental  data. 
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The  aim  of  the  first  part  of  this  paper  is  to  present  an  explicit  analytic  solution  of 
the  residual  stress  problem.  This  analysis  is  then  used  in  conjunction  with  conventional 
linear  elastic  fracture  mechanics  to  predict  the  critical  grain  size. 

It  is  important  here  to  note  a  recent  paper  by  Tvergaard  and  Hutchinson  (1988) 
which  addresses  the  problem  of  residual  stress  microcracking  within  the  framework  of 
thermal  expansion  and  elastic  anisotropy.  For  obvious  reasons,  it  is  difficult  to 
compare  the  predictions  of  the  analysis  presented  here,  which  is  based  on  isotropic 
elastic  response,  with  the  Tvergaard  and  Hutchinson  (1988)  work.  But  some 
discussion  of  the  respective  predictions  is  given. 

Finally,  we  wish  to  draw  attention  to  the  work  of  Fredrich  and  Wong  (1986) 
which  addresses  the  related  problem  of  thermally-induced  microcracking  of  rocks. 
These  authors  assume  elastic  isotropy  together  with  thermal  anisotropy  and  provide 
a  thorough  analysis  of  the  two-dimensional  model  consisting  of  a  square  inclusion 
embedded  in  an  isotropic  aggregate. 

In  this  paper,  we  follow  Evans  (1978)  and  Fu  and  Evans  (1985),  by  considering  a 
plane  hexagonal  grain  array  embedded  in  an  infinite  isotropic  elastic  matrix  (the 
effective  polycrystalline  aggregate).  The  orientation  of  the  axes  of  thermal  anisotropy 
of  the  individual  grains  is  allowed  to  vary.  In  Section  2,  we  give  the  general  solution 
of  the  residual  stress  problem  using  complex  variable  methods.  Amongst  other  things, 
the  analysis  confirms  a  result  due  to  Evans  (1978)  that  the  stress  singularity  at  a  triple 
point  is  logarithmic.  We  go  on  to  discuss  the  stress  intensity  factors  at  putative  cracks 
on  grain  boundary  facets.  It  is  especially  noteworthy  that  we  are  able  to  evaluate  all 
the  integrals  analytically  and  thus  reduce  the  determination  of  stress  intensity  factors 
to  the  evaluation  of  an  algebraic  sum. 

To  determine  the  residual  stress  field  in  an  array  with  given  orientation  distribution 
of  the  axes  of  thermal  anisotropy,  it  is  necessary  to  add  the  contributions  due  to  each 
facet.  The  question  therefore  arises  as  to  how  many  grains  need  be  considered  in  order 
to  get  accurate  results.  It  turns  out  that  it  is  essential  to  consider  at  least  200  grains 
to  get  proper  accuracy. 

Since  the  orientation  distribution  of  the  axes  of  thermal  isotropy  is  not  known,  it 
is  argued  that  the  appropriate  model  for  spontaneous  microcracking  of  randomly 
oriented  arrays  is  to  take  the  most  extreme  local  orientation  at  the  considered  interface, 
together  with  the  ensemble  average  over  all  other  grains.  This  interpretation  shows 
that  one  can  interpret  the  2-grain  and  4-grain  models  of  Fu  and  Evans  (1985)  in  a 
new  light. 

In  Section  4  we  use  the  results  of  our  exact  stress  analysis  to  predict  the  minimum 
facet  size  for  microcracking  during  cooldown.  We  show  that  the  factor  Q  in  (1)  is 
very  sensitive  to  the  assumed  length  of  the  inherent  flaw  in  the  solid.  Finally,  we 
compare  the  predictions  of  this  analysis  with  those  of  other  workers  and  with  some 
experimental  data. 


2.  Analysis 


A  standard  model  of  a  polycrystalline  solid  consists  of  a  regular  hexagonal  array, 
see  Fig.  1 .  In  the  study  of  microcracking  due  to  residual  stresses,  it  is  usual  to  assume 
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Fig.  1.  Hexagonal  grain  array.  Also  shown  are  the  global  and  local  coordinates  Tor  the  rth  grain. 


elastic  isotropy  and  thermal  anisotropy  of  the  individual  grains  with  the  proviso  that 
the  coefficient  of  thermal  expansion  in  the  .x3  direction  of  Fig.  1  is  the  same  for  every 
grain.  At  the  time  of  writing,  little  or  no  data  is  available  concerning  the  orientation 
distribution  of  the  axes  of  thermal  anisotropy.  Thus, we  are  led  to  the  study  of  models 
with  various  orientation  distributions.  However,  throughout  this  paper,  we  will  focus 
on  the  micro-geometry  shown  in  Fig.  1.  We  will  show  that  the  determination  of 
residual  stresses  can  be  accomplished  for  arbitrary  orientation  distributions  of  the  axes 
of  thermal  anisotropy.  Accordingly,  we  defer  the  introduction  of  specific  orientation 
models  to  Section  3. 

The  residual  stress  problem  is,  perhaps,  most  readily  discussed  with  the  help  of  the 
notation  of  Laws  (1973).  Thus,  for  the  rth  grain,  the  constitutive  equation  is 

s'  =  A/«r'+0m',  (2) 

where  s'  and  a'  are  respectively  the  strain  and  stress  in  the  rth  grain.  M  is  the  common 
(isotropic)  compliance  tensor  and  9  the  increase  in  temperature  from  the  stress-free 
configuration.  Also,  m'  is  the  tensor  of  coefficients  of  thermal  expansion. 

When  referred  to  the  local  axes  of  thermal  symmetry,  0.x',.x':.x3,  see  Fig.  1,  the 
tensor  m'  is  expressible  in  the  form 

/a,  0  0\ 

m'  =  0  a2  0  •  (3) 

\0  0  a,/ 

The  components  of  m'  with  respect  to  the  global  axes  O.X|.x,.x,  are  then  found  by 
tensor  transformation : 
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Fig.  2.  Grain  array  showing  those  grains  which  are  considered  in  the  various  /V-grain  models. 


nfu  =  a,  cos2  0,  +a2  sin2  0,, 
m'u  =  (a2— a,)sin  0,  cos  0„ 
nfZ2  =  a,  sin20r+a2cos20r, 

Wu  »  nf2 j  =  0,  =  a3.  (4) 

For  the  polycrystalline  aggregate  as  a  whole,  it  is  obvious  that  the  macroscopic 
compliance  is  M.  Also,  it  is  not  difficult  to  show  that  the  macroscopic  tensor  of 
coefficients  of  thermal  expansion,  m,  is  just  the  orientation  average  of  m'  taken  over 
all  grains : 

m  =  {in'}.  (5) 

With  the  help  of  (4),  it  is  easy  to  calculate  the  overall  coefficient  of  thermal  expansion 
for  arbitrary  orientation  distributions  of  the  axes  of  thermal  anisotropy. 

Throughout  this  paper,  we  will  assume  that  the  polycrystalline  aggregate  has  macro¬ 
scopic  thermal  isotropy  so  that 

mu  =  m22  =  mn  =  a,  (6) 

m23  =  m}[  =  m,;  =  0. 

Nevertheless,  we  note  in  passing  that  it  is  easy  to  extend  the  analysis  to  textured 
aggregates. 

The  residual  stresses  in  the  microcrack-free  solid  maybe  calculated  using  the 
Eshelby  (1961)  technique.  This  calculation  is  performed  by  considering  the  configur¬ 
ation  of  Fig.  1  (see  Fu  and  Evans,  1985)  wherein  the  hexagonal  array  is  surrounded  by 
a  material  whose  properties  are  those  of  the  effective  aggregate.  The  number  of  grains 
considered  is  at  our  disposal.  Thus,  for  given  orientation  of  the  thermal  axes  of  an 
array,  we  can  have  a  hierarchy  of  models.  To  emphasize  the  point,  let  us  suppose  that 
we  are  interested  in  the  interfacial  stresses  on  a  particular  facet,  AB.  in  Fig.  2.  In  the 
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first  instance,  we  can  consider  the  2-grain  model  consisting  of  the  two  grains  labelled 
1  surrounded  by  the  effective  aggregate.  Next,  we  can  consider  a  4-grain  model 
consisting  of  the  two  grains  labelled  1  plus  the  two  grains  labelled  2  surrounded  by 
the  effective  aggregate.  As  shown  in  Fig.  4,  we  can  then  consider  the  2+2+6  =  10- 
grain  model,  followed  by  the  2+2+6+12  =  22-grain  model,  etc.  Clearly,  other 
models  for  the  given  orientation  are  possible.  However,  as  the  considered  number 
of  grains  increases,  the  predicted  stresses  on  the  facet  AB  will  tend  to  a  limit.  An 
important  part  of  this  investigation  is  to  determine  the  number  of  grains  which  must 
be  considered  in  order  to  get  accurate  estimates  of  the  stresses  on  the  facet  AB. 

Returning  now  to  the  general  problem,  we  recall  that  the  Eshelby  method  requires 
that  each  grain  and  the  surrounding  matrix  be  considered  as  isolated  and  the  tem¬ 
perature  increased  by  9.  Thus,  the  rth  grain  is  now  subject  to  strain  9mr.  We  next 
strain  each  grain  so  that  it  has  the  same  strain  as  the  surrounding  matrix  and  all  grains 
can  then  fit  together  into  the  matrix.  The  required  extra  strain  of  the  rth  grain  is 

el  =  fl(m— m')  (7) 

and  this  is  accomplished  by  the  application  of  surface  stresses  s£  where 

s;  =  U,  (L  —  M~').  (8) 

The  final  step  in  the  procedure  is  to  superpose  a  layer  of  body  force  over  the  surface 
of  each  grain  to  annihilate  the  surface  tractions  obtained  from  (8).  Thus,  on  the 
interface  between  the  rth  grain  and  the  <?th  grain,  the  resultant  superposed  body  force 
is 


t=  -sjn— s?(  — n) 
*«-<)■ 


=  0L(mr-m,)n,  (9) 

where  n  is  the  unit  normal  from  the  rth  grain  to  the  <yth  grain.  We  note  that  the 
traction  vector  t  is  constant  over  the  interface.  It  then  follows  that  the  entire  residual 
str  .  -s  distribution  in  the  uncracked  body  is  obtained  by  superposition  of  the  stresses 
<  a  ’the  stresses  induced  by  the  totality  of  interfacial  body  forces  (9).  It  is  important 
to  rei.iember  the  contribution  to  the  stress  field  due  to  the  interface  between  the  outer 
grains  and  the  surrounding  matrix. 

It  is  therefore  clear  that  the  residual  stress  problem  for  the  microcrack-free  aggregate 
is  solved  once  we  can  evaluate  the  stress  distribution  in  an  infinite  solid  due  to  constant 
body  force  over  a  single  line  segment — superporition  gives  the  final  solution.  This 
problem  is  readily  solved  using  complex  variable  methods. 

First,  consider  a  point  (or,  more  precisely,  line)  force  with  components  P ,,  Pz  at 
the  point  represented  by  the  complex  number  :0  =  .x0  +  iy0  in  an  infinite  isotropic 
solid.  For  this  problem,  it  is  well  known  that  the  complex  potentials  are  given  by 


0(-)  =  - 


2ji(k+  1) 


ln(r---0). 


*(->  = 


kP 


2tc(k  +  l) 


In  (:-z0)+z0 


2ji(k  +  l) 


(10) 

(ID 


608 


N.  Laws  and  J.  C.  Lee 


where  P  =  P,  +iP2.  Also,  since  we  are  interested  in  plane  strain,  k  =  3-4v,  but  the 
work  described  below  also  applies  to  generalized  plane  stress  when  k  =  ( 3  —  v)/(  1  +  v). 
Next,  consider  the  straight  line  segment  joining  two  vertices  (i.e.  triple  points)  r, 
and  r:,  and  suppose  that  this  segment  (grain  boundary  facet)  is  subject  to  uniform 
distributed  body  force  t  =  t ,  +  i/2  per  unit  length.  It  may  be  shown  that  differentiation 
of  (10)  and  (11)  with  respect  to  r,  followed  by  integration  along  the  segment  r 
implies  that 


<!>(-)  =  <t>'(z) 


27t(»C+  l)co 


(12) 


V(Z)  =  *'(=)  = 


0cfa>+/c3)|c£>|  z— r2 
27t(K+l)tU2  0Z  —  2| 


fM  Zi 
27t(K+  1)W  L-  — -2 


(13) 


where  as  -  r2  — r,.  It  is  now  possible  to  obtain  the  complex  potentials  for  the  entire 
distribution  of  body  forces  over  grain  boundary  facets  and  thus,  the  residual  stress 
field  in  the  uncracked  solid.  While  the  algebraic  details  are  too  cumbersome  to  be  fully 
reported  here,  some  particular  considerations,  related  to  the  stress  singularities  at  the 
triple  points,  are  useful. 

From  (12)  and  (13)  the  total  complex  potentials  are  clearly  of  the  form 

$(2)-  £  Ak  In  (z—zk),  (14) 

k-  I 


*v  Akf, 


k+k 


¥(2)=  £  5*  In  (2-2*)-  £ 

*•1  k-\  +  ~  -* 


(15) 


where  zk(k  =  1,2,..  .N)  is  the  set  of  triple  points  and  where  Ak,  Bk  are  a  set  of  known 
constants.  With  the  help  of  the  standard  formulae 

Vxx  +<ry,  =  2{4>(r)+«I>(r)}, 

<T,r  +  iff,,  =  0(2)  +  +  2<I>'  (2)  +  ¥(2), 


we  can  compute  the  stresses  in  the  form 


ff  XX  4*  ff  „ 


=  2  [£  [Ak  In  (z  -zk)  +  Ak  In  (f  -  i*)] j ,  (16) 


=  1  [Ak  +  Bk)  In  (2 — zk) + Ak  In  (f — 2*)] 4-  £  -fll ,  (i7) 

*  -  I  k-\  -~-k 


Equation  (17)  immediately  suggests  the  possibility  of  an  r" 1  type  singularity  at  each 
triple  point.  However,  if  we  write 


2-2*  =  r*  exp  (i#*), 


then  it  is  clear  that 
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Fig.  3.  Grain  boundary  crack  of  length  2 a  at  the  distance  b  from  the  triple  point  A. 


7- r  =  exp(-2i0*). 

Hence  there  is  no  r” 1  type  singularity  and  the  only  stress  singularity  is  logarithmic. 
Obviously,  this  singularity  exists  at  every  triple  point,  provided  the  axes  of  thermal 
anisotropy  of  the  three  associated  grains  do  not  coincide. 

It  is  abundantly  clear  from  (9),  (16)  and  (17)  that  evaluation  of  the  global 
residual  stress  field  has  been  reduced  to  an  algebraic  problem.  However,  the  complexity 
of  the  algebra  should  not  be  underestimated. 

Let  us  now  turn  to  the  evaluation  of  stress  intensity  factors  which  are  relevant  in 
the  study  of  grain  boundary  microcracking  due  to  residual  stress.  It  is  commonly  held 
that  grain  boundary  microcracking  is  due  to  triple  point  defects.  But  there  remains 
the  possibility  that  defects  occur  elsewhere  on  the  grain  boundary.  Since  the  stress 
intensity  factors  for  arbitrarily  located  grain  boundary  cracks  are  readily  inferred 
from  our  previous  analysis,  we  consider  the  putative  crack  CD  on  the  facet  AB  as 
shown  in  Fig.  3.  For  the  required  calculations  it  is  convenient  to  orient  the  global 
x,-axis  along  that  facet.  It  then  follows  that  the  interfacial  normal  and  shear  stresses 
are  given  directly  by  (17). 

A  standard  result  in  linear  elastic  fracture  mechanics  shows  that  the  stress  intensity 
factors  at  D  are  obtained  from 


K{  +  ii^H  =  — I  (av,+iaxy)  - -dr.  (18) 

V  a  1 

From  (17)  and  (18)  we  observe  that  the  calculation  of  Kt  and  Kn  involves  the 
evaluation  of  some  elementary  integrals  together  with  integrals  of  two  types : 


J . 
J2 


i 

f 


fl+t  1 
1-r  r+c 


dr. 


r  j—  In  (r+c)  dr. 


where  c  is  complex  and  |c|  ^  1.  The  values  of  7,  and  7:  may  be  obtained  in  closed 
form  as 


Fig.  4.  Grain  array  with  regular  orientation  distribution. 


With  the  help  of  these  explicit  expressions  for  /,  and  /j,  it  follows  that  the  deter¬ 
mination  of  stress  intensity  factors  in  the  aggregate  has  been  reduced  to  a  simple 
algebraic  sum.  It  is  neither  illuminating  nor  necessary  to  include  details  of  the  tech¬ 
niques  which  we  use  to  handle  the  orientation  distributions  of  the  axes  of  thermal 
anisotropy  and  the  summation  over  various  grain  arrays. 

We  emphasize  that  the  preceding  discussion  also  shows  that  we  can  readily  compute 
residual  stresses  and  the  associated  stress  intensity  factors  for  grain  boundary  micro- 
cracks  for  an  arbitrary  crystalline  array,  regular  or  otherwise. 


3.  The  Residual  Stress  Field 

The  analysis  presented  in  preceding  sections  provides  an  analytic  solution  for  the 
entire  residual  stress  field  in  the  polycrystalline  aggregate.  From  a  practical  point  of 
view  the  usual  thermoelastic  moduli  of  an  individual  grain  are  taken  as  data,  but  the 
orientation  distribution  of  the  various  grain  arrays  is  not  available  at  the  time  of 
writing.  Thus,  it  will  be  essential  to  quantify  the  effect  of  grain  orientation  on  the 
residual  stress  field  and  hence  on  microcracking  during  cooldown.  However,  from  a 
computational  standpoint,  the  first  issue  to  be  addressed  must  be  the  determination 
of  the  number  of  grains  which  is  needed  to  get  an  accurate  value  for  the  interfacial 
stresses.  _ 

For  illustration,  we  consider  the  regular  orientation  model  shown  in  Fig.  4.  First 
we  remark  that  for  this  model,  the  overall  coefficients  of  thermal  expansion  may  be 
obtained  from  (5)  as 


mu  =  m22  *  («i  +ctj)/2,  m}}  =  a3. 

At  this  stage,  it  is  helpful  to  introduce  some  notation  which  is  consistant  with  that  of 
Evans  (1978).  Thus,  let  the  temperature  drop  during  cooldown  be  AT,  let  the  differ¬ 
ence  between  the  maximum  thermal  expansion  coefficient  and  the  average  be  denoted 
by  Act.  and  let  /  be  the  facet  length. 
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Fig.  5.  Residual  stress  at  the  midpoint  of  AB  vs  number  of  grains.  Note  the  variable  scale  on  the  .t-axis. 


We  successively  compute  the  normal  and  shear  stresses,  on  AB,  for  the  2-grain 
model,  the  4-grain  model,  etc.  The  predicted  value  of  the  normal  stress  at  the  midpoint 
of  AB,  for  the  various  multi-grain  models,  is  shown  in  Fig.  5.  We  note  that  it  is 
essential  to  consider  an  array  consisting  of  at  least  200  grains  to  get  an  accurate 
estimate  of  the  normal  stress  at  the  midpoint  of  the  facet.  We  cannot  emphasize 
too  strongly  that  the  preceding  assertion  also  applies  to  other  fixed  orientation 
distributions.  For  example,  we  have  investigated  the  grain  orientation  distribution 
studied  by  Tvergaard  and  Hutchinson  (1988),  and  we  have  found  that  the  results 
obtained  using  our  technique  with  268  grains  are  entirely  consistent  with  the  finite 
element  results  of  these  authors.  Thus,  we  repeat,  that  in  order  to  get  accurate 
predictions  of  the  residual  stresses  on  a  particular  facet  in  a  given  polycrystal,  it  is 
necessary  to  consider  at  least  200  grains  surrounding  that  facet. 

It  is  of  interest  to  use  Fig.  6  to  calculate  the  residual  normal  stress  on  the  facet  AB 
for  the  particular  alumina' AD995.  Here  we  can  take 


Fig.  6.  Residual  stress  along  the  grain  boundary  AB  for  various  grain  arrays. 
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Fig.  7.  Geometries  and  orientations  for  (a)  2-grain  model  and  (b)  4-grain  model. 


E  -  372  GPa,  v  =  0.25, 

Aa  =  10-6C~ ',  AT=  1500°C, 

and  it  follows  that  the  normal  stress  at  the  midpoint  of  AB  is  150  MPa.  This  result 
may  be  contrasted  with  the  tensile  strength  for  AD995  which  is  262  MPa. 

An  equally  important  issue  pertains  to  sensitivity  of  the  interfacial  stresses  to  the 
assumed  orientation  distribution  of  the  axes  of  thermal  anisotropy.  As  has  been 
emphasized,  the  orientation  distribution  function  for  a  typical  polycrystal  is  not 
known.  Indeed,  it  is  evident  that  the  orientation  distribution  will  vary  from  sample 
to  sample.  Therefore,  when  we  examine  models  in  which  the  orientation  distribution 
is  not  fixed  a  priori,  it  is  clear  that  we  need  to  consider  the  appropriate  ensemble 
average.  In  this  sense  we  can  arrive  at  an  alternative  interpretation  of  the  various 
multi-grain  models. 

To  explain,  let  us  focus  on  the  facet  which  will  experience  the  greatest  residual 
normal  stress.  We  therefore  consider  two  grains  whose  maximum  contraction  will  be 
perpendicular  to  their  interface,  see  Fig.  7(a).  Now  regard  these  grains  as  fixed.  Since 
the  orientation  of  the  remaining  grains  is  random,  we  can  compute  the  ensemble 
average  of  the  interfacial  stresses  by  considering  a  large  number  of  samples  with 
randomly  oriented  axes  of  thermal  anisotropy.  But  it  is  evident  that  the  model  obtained 
by  taking  two  fixed  grains  together  with  the  ensemble  average  over  all  orientations  of 
the  surrounding  grains  is  equivalent  to  the  model  in  which  the  two  fixed  grains  are 
surrounded  by  the  effective  aggregate.  We  refer  to  this  model  as  the  2-grain  model 
for  the  randomly  oriented  aggregate.  This  model  will  give  us  the  average  stress  on  the 
interface  AB.  In  similar  fashion  we  can  consider  the  model  associated  with  the  4  fixed 
grains  of  Fig.  7(b).  The  model  in  which  these  4  grains  are  surrounded  by  the  effective 
aggregate,  will  give  us  the  average  stress  on  AB — the  average  referring  to  all  orien¬ 
tations  of  the  remaining  grains.  In  this  way  we  can  arrive  at  the  W-grain  model  for 
the  randomly  oriented  aggregate. 

It  is  obvious  that  the  matheihatics  of  the  iV-grain  model  for  the  aggregate  with 
given  a  priori  orientation  distribution  is  the  same  as  the  mathematics  of  the  iV-grain 


613 


Microcracking  in  polycrysutUine  ceramics 


tmi  ( 

Fig.  8.  Mode  I  stress  intensity  factors  for  cracks  variously  located  in  the  facet  AB :  regular  array  orientation 

model. 

model  for  the  randomly  oriented  aggregate.  However,  the  physical  interpretation  of 
the  two  models  is  completely  different.  To  elaborate,  we  refer  to  Fig.  6,  which  amongst 
other  things,  shows  the  predicted  normal  stress  for  the  2-grain  model  of  the  given 
regular  orientation  distribution  of  Fig.  4.  But  we  can  also  interpret  the  2-grain  model 
curve  in  Fig.  6  as  giving  the  average  normal  stress  on  AB  for  a  polycrystalline 
aggregate  wherein  two  grains  are  held  fixed  and  the  remainder  allowed  to  take  all 
orientations. 

While  the  magnitudes  of  the  residual  interfacial  stresses  are  of  interest  and  sig¬ 
nificance  in  their  own  right,  it  is  perhaps  more  relevant  to  consider,  in  greater  depth, 
the  stress  intensity  factors  which  would  exist  at  the  tips  of  interfacial  cracks.  Indeed 
a  careful  analysis  of  the  various  stress  intensity  factors  is  an  essential  prerequisite  to 
any  understanding  of  microcracking  due  to  cooldown  or  subsequent  mechanical 
loading. 

We  first  consider  the  regular  array  orientation  of  Fig.  4  and  examine  the  stress 
intensity  factors  for  variously  located  cracks  on  the  interface  AB,  as  shown  in  Fig.  3. 
We  emphasize  that  for  this  configuration,  the  directions  of  maximum  contraction  of 
the  grains  on  either  side  of  AB  are  normal  to  the  interface.  The  mode  I  stress  intensity 
factor  is  shown  in  Fig.  8  as  a  function  of  initial  defect  length  for  various  locations  of 
the  defect.  We  here  restrict  our  attention  to  2a//  <  0.3.  We  see  that  for  given  initial 
defect  length,  the  most  critical  crack  emanates  from  the  triple  point.  In  addition,  we 
emphasize  that  Kf  is  sensitive  to  the  length  of  the  initial  defect. 

Figures  9  and  10  show  the  variation  in  stress  intensity  factors  at  the  right  hand  tip 
for  triple  point  cracks  on  AB  within  the  regular  orientation  array  of  Fig.  4.  except 
that  the  orientation  of  one  grain  is  allowed  to  vary.  We  note  the  sensitivity  of  the 
results  to  orientation  of  the  single  grain  at  the  left  of  the  critical  triple  point.  In  the 
same  spirit.  Figs  1 1  and  12  show  the  effect  of  changing  the  orientation  of  one  grain 
to  the  right  of  B  in  the  otherwise  regular  array  of  Fig.  4.  In  either  case.  Kn  is  much  less 
than  K,.  Finally,  we  investigate  the  effect  of  changing  the  orientation  of  the  single 
grain  above  AB.  The  effect  of  orientation  on  and  Ku  as  a  function  of  orientation . 
is  shown  in  Figs  1 3  and  1 4.  We  note  that,  even  for  small  defects,  the  effect  of  orientation 
can  be  significant. 
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Fig.  9.  Mode  I  stress  intensity  factors  for  various  orientations  of  the  grain  to  the  left  of  the  triple  point 
crack — in  the  otherwise  regular  orientation  anay. 
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Fig.  10.  Mode  II  stress  intensity  factors  for  various  orientations  of  the  grain  to  the  left  of  the  triple  point 
crack — in  the  otherwise  regular  orientation  array. 
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Fig.  1 1.  Mode  I  stress  intensity  factors  for  various  orientations  of  the  grain  to  the  right  of  the  triple  point 
crack — in  the  otherwise  regular  orientation  array. 
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Fig.  12.  Mode  II  stress  intensity  factors  for  various  orientations  of  the  grain  to  the  right  of  the  triple  point 
crack — in  the  otherwise  regular  orientation  array. 


Fig.  13.  Mode  I  stress  intensity  factors  for  various  orientations  of  the  grain  above  the  triple  point  crack — 

in  the  otherwise  regular  orientation  array. 


Fig.  14.  Mode  II  stress  intensity  factors  for  various  orientations  of  the  grain  above  the  triple  point  crack — 

in  the  otherwise  regular  orientation  array. 
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Fig.  IS.  Mode  I  stress  intensity  factors,  for  triple  point  cracks,  for  various  grain  arrays. 


It  is  inappropriate  to  produce  a  multitude  of  figures  to  show  the  sensitivity  of  the 
stress  intensity  factors  to  orientation  changes.  Nevertheless,  we  can  assert  that,  for 
the  regular  orientation  array,  microcracking  is  essentially  a  mode  I  phenomenon  and 
that  the  orientation  of  contiguous  grains  is  important. 

Next,  we  consider  some  problems  associated  with  the  possibility  of  random  orien¬ 
tation  of  the  axes  of  thermal  anisotropy.  Here  we  focus  on  the  2-grain  and  4-grain 
models  in  Fig.  7  which  were  discussed  earlier.  Hence,  these  models  will  give  us  ensemble 
averages  for  the  stress  intensity  factor.  Figure  15  shows  the  stress  intensity  factors 
which  are  predicted  by  the  two  models.  For  comparison,  the  regular  orientation  array 
result  is  also  included  in  Fig.  15.  Of  course,  for  the  regular  array  Kn  =  0. 

The  obvious,  and  disconcerting,  conclusion  is  that  the  various  models  give  sig¬ 
nificantly  different  results.  The  status  of  the  regular  orientation  array  is  clear.  And  it 
is  likely  that  the  4-grain  model  of  Fig.  7(b)  gives  the  most  extreme  local  orientations 
for  potential  cracking  on  the  interface  AB.  Thus,  when  we  consider  grain  boundary 
microcracking,  we  compare  the  regular  orientation  model  with  the  4-grain  model. 
The  reader  is  warned  that  other  models  give  different  results. 


4.  Grain  Boundary  Microcracking 

The  prediction  of  initial  microcracking  during  cooldown  is  achieved  by  focusing 
on  the  worst  possible  scenario.  Thus,  we  consider  triple  point  cracks  within  the  regular 
orientation  array  or  the  4-grain  model  of  Fig.  7(b).  We  recall  that  for  practical 
purposes,  AT,,  is  negligible.  From  Fig.  15  we  have 

^W(^).  cm 

pAaATy/l  V  1  > 

where  /depends  upon  the  considered  model.  Next  let  G|(l  be  the  grain  boundary, 
mode  I,  plane  strain  toughness.  It  then  follows  from  (19)  that  microcracking  during 
cooldown  first  occurs  at  critical  grain  size,  l„  where 
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(20) 


where  2a  is  the  length  of  the  inherent  triple  point  crack.  An  important  parameter  here 
is  the  ratio  (2aile).  As  is  obvious  from  Fig.  15,  /  changes  rapidly  for  0  <  2a/ 1  <0.1. 
It  is  therefore  clear  that  any  prediction  of  critical  grain  size  is  highly  sensitive  to  the 
flaw  size  in  the  material. 

Let  us  now  discuss  some  critical  grain  size  predictions  for  specific  polycrystalline 
ceramics.  In  the  first  instance  we  consider  the  magnesium  titanate  system  studied  by 
Evans  (1978).  For  this  system  Evans  (1978)  suggests  2a//  =  0.1.  We  can  then  read 
off  the  value  of /(0.1)  from  Fig.  15 : 

y (0.1)  *  2.85  (regular  orientation  array) 

=*  3.99  (4-grain  model). 

If  we  take  the  value  v  =  0.25,  the  above  formulae  in  turn  lead  to 


4 


=  7.39 

=  4.02 


E(A*AT)2 

E(AaAT)2 


(regular  orientation  array). 


(4-grain  model). 


The  most  recent  formula  proposed  by  Fu  and  Evans  (1985)  is 

1(1 +  v  )2Gtb 
k  E(AadT)2  ' 


(21) 


The  factor  /?  is  chosen  to  be  approximately  3.5  in  order  to  get  an  adequate  correlation 
with  experiment.  If  we  interpret  correctly,  the  proposed  formula  is  independent  of  the 
ratio  (2a//).  For  v  =  0.25  the  Fu  and  Evans  (1985)  formula  becomes 


4 


5  47 - - 

E(AaAT)2  ’ 


To  calculate  the  critical  grain  size,  we  use  the  data  given  by  Evans  (1978),  namely 

£  =  240  GPa,  Act  =  5  x  IO^C-1,  AT  =  1000°C. 

In  addition,  Evans  (1978)  suggests  that  the  value  of  Gtb  lie  between  2  Jm~:  and  10 
Jm”  2.  The  various  predictions  of  the  three  models  are 


4  =  4.9-24.6  gm  (regular  orientation  array) 

=  2.7-1 3.4  ^m  (4-grain  model) 

=  3.6-l8.3/im(Fu  and  Evans). 

The  critical  average  grain  size  inferred  from  experiment  [2]  is  about  3  gm,  but  as 
emphasized  by  Evans  (1978),  the  experimental  value  must  be  interpreted  with  care. 
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In  the  same  spirit,  we  wish  to  emphasize  the  critical  role  played  by  the  assumption 
that  (2a/ 1)  =  0.1. 

Next  we  turn  to  the  A1,03  studied  by  Rice  andPoHANKA  (1979)  and  by  Tvergaard 
and  Hutchinson  (1988).  Following  these  authors  we  here  take 

E  =  350  GPa,  v  =  0.2,  Aa  =  0.55  x  10" ‘C"1,  AT  =  1000°C, 

2a// =  0.02,  (?|b  =  2Jm-2. 

From  Fig.  14 

/( 0.02)  =  2.21  (regular  orientation  array) 

=  2.63  (4-grain  model). 

It  now  follows  that 

4  =  240  pm  (regular  orientation  array) 

=  1 70  pm  (4-grain  model) 

=  103  ^m(Fuand  Evans). 

Tvergaard  and  Hutchinson  (1988)  also  studied  this  Al,Oj  system  with  the  help 
of  their  periodic  array  and  found  4  =  216  pm.  For  the  same  array,  the  present  model 
gives  4  =  218  pm.  The  experimental  value  is  /c  =  200  pm. 
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The  effect  of  microcracks  on  energy  density 


13-1  Introduction 

la  ci"grit  years,  there  has  been  considerable  Interest  in  the  response  of 
polycrystailine  brittle  solids  which  contain  a  homogeneously  distributed 
family  of  microcracks.  As  far  as  the  author  is  aware,  the  quantification  of  the 
effect  of  such  a  population  of  aricrocracla  on  the  macroscopic  response  of  the 
solid  waa  initiated  by  Budiansicy  and  O'Comxdl  [I].  Other  contributions  have 
been  gtven  in  a  series  of  papers  by  Evans  and  cowodoas  C2-6],  Hoaglasd  and 
Embury  (71,  Clarke  [8],  ffocnig  [9],  Laws  and  Brocken brough  [10,  11] 
Hutchinson  [12]  and  Ctaraiamhides  and  McMeririsg  [13]. 

The  topics  of  major  concent  in  this  literature  are  the  loss  of  stiffness  due 
to  a  given  family  of  microcracks  and  the  audeariou  of  Slither  microcracks. 
In  addition,  there  is  considerable  interest  in  the  extent  to  which  microcracks 
in  the  process  zone  of  a  macroscopic  crack  can  have  a  shielding  effect  an 
the  macroscopic  crack  tip.  These  issues  are  further  discussed  in  this  paper. 
The  development  is  presented  in  the  same  spirit  as  those  of  earlier  studies 
[1-13]. 

All  the  foregoing  contributions  refer  to  open  microcracks.  And,  for  the  most 
part,  the  respective  authors  focus  on  loss  of  macroscopic  stiffness,  rather  than 
change  of  macroscopic  energy  density  -  sines  the  latter  (when  relevant)  is 
trivially  obtained  from  the  former. 

At  this  stage  of  our  discussion  it  is,  perhaps,  essential  to  emphasize  a  point 
which  is  developed  by  Hutchinson  [12]  and  by  Charalambidcs  and 
McMeeking  [13],  namely  that  existence  of  the  standard  microscopic  energy 
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Chapter  13 


density  does  not  imply  existence  of  a.  macroscopic  energy  density  in 
t  tnicgocracking'solid.  In  fact,  as  will  be  shown  lazes,  severe  restrictions  must 
be  placed  on  the  microcode  audeation  function  in  order  that  the  microcode- 
ing  polycrystalline  aggregate  con  be  regarded  as  masoscopicaily  hyperdasac 
In  addition  we  draw  attention  to  the  important  paper  of  Horn  and  Nemat- 
Nasscr  £14]  which  pertains  to  the  etTecs  of  microcode  closure  and  friction  on 
the  macroscopic  response.  These  authors  show  that  the  macroscopic  com¬ 
pliance  tensor  of  such  a  solid  depends  on  the  load  path  and  need  not  even  be 
symmetric 

In  this  paper,  a  summary  of  results  is  given  Gar  loss  of  stiffness  of  a  solid 
containing  various  families  of  microcodes.  For  simplicity,  we  restrict 
attention  to  dilute  dispersions  of  microcodes  (so  that  the  microcracks  do  not 
interact).  Further,  we  make  the  standard  assumption  that  at  least  in  a  first 
approximation,  each  grain  of  the  poly  crystalline  aggregate  is  isotropic  and 
that  residual  stresses  may  be  neglected.  Next,  we.  discuss  some  issues  which 
Gram  the  choice  of  the  microcode  audeation  function.  In  the  absence  of 
definitive  experimental  evidence,  a  simple  audeation  function  is  proposed 
which  leads  to  the  tnicrocracked  solid  being  hypcrciastic  This  choice  of 
audeation  function  is  not  essential  but  Leads  easily  to  an  aiwnimeut  of  the 
urtfmr  of  microcadc  <hiei<iinj  on  a  macroscopic  crack  tip.  The  results  thus 


obtained  are  i 


and  McMeeking  £13],  and  Ortiz  £16]. 


13.2  Mica  cracked  solid  with  given  crack  density 


The  approach  adopted  hero  doseiy  follows  that  ofBudiansky  and  O'Connell 
£VL  Law*  and  Brockenbraugh  £10],  and  Hutchinson  £12].  Consider  an 
isotropic  solid  with  compliance  tensor  AfQ  which  contains  a  family  of  variously 
oriented  open  microcracks.  It  is  essential  to  assume  that  the  micro  cades  are 
homogeneously  dispersed.  And  it  is  convenient  to  assume  that  the  mioocracks 
axe  of  similar  elliptic  piamorm  and  that  crack  orientation  and  size  ate  not 
correlated.  A  convenient  measure  of  micro  crack  density  has  been  suggested  by 
Budiansky  and  O'Connell  £1]: 


where  *Y  is  the  number  of  cades  per  unit  volume.  A  is  the  area  of  a  crack  and 
P  is  perimeter.  Also  (  }  denote  the  orientatioa  average  of  the  bracketed 
quantity. 

From  £10],  we  can  immediately  read  ofT  the  formula  for  the  loss  of  stiffness 
of  a  solid  containing  a  dilute  distribution  of  «mflar  microcracks; 
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Af-M0+JJ^!{A}. 


(113 


where  A  is  a  tensor  whose  components  are  given  in  CIO]  tor  &  variety  of 
microcxack  «hap,^li  and  where  £(£)  is  the  elliptic  integral  of  the 

irrrm'i  kind. 

It  is  of  interest  to  give  some  explicit  deductions  from  equation  (113  (hr 
ptt+icMiar  microcrack  distributions.  Far  a  three-dimensional  randomly 
family  of  penny-shaped  we  can  recover  the  Budiansky  and 

O’Connell  Cl]  results 


£  t  -16(1 -*SttO -3*0) 
£„“  *  45(2- vj 


(133 


0  „  32(1  —  Vq)  (5  —  Vq) 

G*ml~P  45(2  —  Vq)  ’ 


(13.4) 


where  G  is  the  shear  modulus  and  v  is  Poisson’s  ratio.  Also,  lor  a  family  of 

gfifliw-rf  pj»tmy-«hap»H  cracks 


|.i->-*3 


03-5) 


C  .  16al-*e 


where  £  is  now  Young’s  modulus  normal  to  the  microcrack  &ces  and  Gis  the 
transverse  shear  modulus. 

When  we  have  a  family  of  aligned  slit  cracks,  we  find  that  the  reduction  in 
Young’s  modulus  normal  to  the  slits  is  given  by 


(1 3.7) 


whereas  for  a  two-dimcnsionally  randomly  oriented  family  of  slits  [10] 
£o  4 

Other  explicit  results  are  obtainable  but  details  are  not  given  here. 


(133 
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L3-3  Micro  crack  audeadoe 

Some  potential  audeadon  functions  have  been  discussed  by  Brockenbrough 

Ortiz  £16].  It  is  particularly  appropriate  here  to  consider  the  audeadon 
fancdonin£l2]  which  is  intimately  associated  with  an  isotropic  distribution  of 
tfaree-dimeasionai  randomly  oriented  penny-shaped  crarira  According  to 
Cb+niamkiAeai  and  McMeeicing  £13],  a  good  approximation  to  the  Budian- 
sky  and  O’Connell  £1]  estimates  for  trifThna  loss  is  given  by 


«  »  . _.i£ 

l£m7.mt  *7- 


(13-3) 


far  all  crack  less  than  9/16.  The  teasorial  stress-strain  relation 

corresponding  to  equations  of  (119)  is  given  by  Hutchinson  £12]  in  the  form 


(mo) 


Further,  the  associated  audeadon  fancrion  is  proposed  Gar  proportional 
loading  to  be 


(13.11) 


where 


p-ia-a)1 


(13-12) 


Thus,  under  proportional  loading,  0  is  taken  to  increase  monotonically  until  it 
reaches  a  saturation  value  £r  It  then  follows  from  equations  (13-10),  (13.11) 
and  (13.12)  that  a  macroscopic  energy  density,  U,  exists.  In  fact, for  the  dilure 


CT«|r-M0<r  +  ~F(p) 
where 

F{p)  -  gf[p\ 
so  that 

3U 

*»r-. 


(13.13) 

(13.14) 

(13.15) 
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While  it  is  true  that  equation  (1X9)  gives  a  good  approximation  to  the 
Budiansky  and  O'Cooaril  [1}  results  over  a  large  range  of  crack  densities,  it  is 
equally  true  that  the  approximations  embodied  in  equation  (1X9)  are  sot 
necessarily  good  for  small  values  of  U  -  as  is  evident  from  equations  (13  J)  and 
(1X4).  This  observation  has  prompted  the  investigation  of  a  different 
nudearion  function  which  dosdy  approximates  the  form  in  equation  (13X1) 
but  which  also  leads  to  of  a  macroscopic  energy  density  function. 

The  nudearion  Auction  proposed  for  investigation  here  is  obtained  by 
demanding  that  the  macroscopic  stress-strain  relation  in  equation  (1X2)  be 
derivable  Dorn  a  macroscopic  energy  density  U(a).  Indeed,  it  is  not  difficult  to 
show  thin  requirement  ^mpii***  the  nudeatjon  Auction 
proportional  loading  is 

-  9(t\  (1X16) 

where 

«*-*•{ A}<r.  (1X17) 

Also  the  macroscopic  energy  density  is  here 
l  g£fJc) 

0'«±<r-M0<r+-^S<Ht).  (1X18) 

where 

<7(1)  -  cy(r).  (1119) 

dearly  the  nudearion  Auction  in  equation  (1X16) « inthnateiy  bound  up  mith 
the  relation  in  equation  (13.2)  vice- versa.  At  juncture,  we 

are  content  to  observe  that  the  proposed  fanewVyr  in  equation  (1X16)  g*hfhrf« 
the  same  anisotropy  as  the  stress-strain  relation. 

We  note  that  other  nudearion  Aucsiops  which  do  got  lead  to  a  macroscopic 
energy  density  Auction  have  been  considered  by  Hutchinson  £12]. 

A  particularly  useful  consequence  of  the  of  an  energy  fundion  is 

that  we  can  easily  determine  the  of  shielding  on  a  stationary 

macroscopic  crack  due  to  microcracking  in  the  process  zone.  This  is  readily 
accomplished  using  the  J  integral.  Following  Hutchinson  [12],  we  only 
consider  Mode  I  small-scale  n^pgaring  an  argument  ofEvans 

and  Faber  [J],  Laws  and  Brocken  brough  [11),  and  Hutchinson  [12],  we  assert 
that  Aar  a  contour  remote  from  a  macroscopic  crack  in  an  isotropic  solid 

(l-rSlC1 

£0  * 


J 


0X20) 
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where  1C  is  the  "applied”  stress  intensity  factor. 

On  the  other  hand,  for  contours  in  the  saturation  none, 

,  a -»***«. 


(U21) 


Invariance  of  J  gives  a  formula  for  for  isotropic  materials. 

Thus  for  a  randomly  oriented  distribution  of  pemxy-sfcaped  cracks  we  Sad 


¥-1-«fe10-3’«+»’5-3* 


(1122) 


For  anisotropic  distributions  of  microcracks,  the  algebraic  details  are  more 
tedious.  Omitting  details,  we  assert  that  for  a  distribution  of  aligned 
pennyshaped  cracks  in  the  process  none 


£hr  t  1<g-7vM 
k  “fia-voT** 

Also  for  aligned  slit  stacks 
K  "l  T? 

whereas  for  two-dimensinnafly  randomly  oriented  slit  cracks 

£*ml 
k  1  a  *** 


(1123) 


(1124) 


(1325) 


The  preceding  formula  axe  readily  compared  and  contrasted  with  the  formulae 
obtained  for  different  nudcarion  ftacrions  by  Hutchinson  [12],  Chazaiam- 
bides  and  McMeddng  [12]  and  Ortiz  [16]. 
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ABSTRACT 

In  this  paper  we  consider  the  problem  of  a  debond  crack  on  the  interface  of  a  circular  fiber 
in  an  infinite  matrix,  when  the  matrix  is  loaded  by  uniaxial  tension  at  infinity.  We  pay  special 
attention  to  those  orientadons  of  applied  load  for  which  the  crack  remains  open.  Finally  we  use  the 
analysis  to  calculate  the  loss  of  stiffness  in  an  unidirectional  ceramic  matrix  composite  due  to 
interfacial  microcracfcing. 

INTRODUCTION 

This  paper  is  concerned  with  fiber-matrix  interfacial  debonds  in  a  unidirectional  fiber 
reinforced  composite.  In  the  first  place  we  discuss  the  problem  of  a  single  debond  crack  at  the 
interface  between  a  circular  fiber  and  an  infinite  matrix.  The  crack  is  assumed  to  be  open  with  zero 
tractions  on  the  crack  faces  and  the  matrix  is  loaded  by  uniaxial  tension  at  infinity.  This  problem 
has  already  been  discussed  by  England  (1966),  Perlman  and  Sih  (1967),  Toya  (1974)  and  Piva 
(1982).  We  draw  explicit  attention  to  the  fact  that  the  crack  closes  for  a  significant  range  of 
orientations  of  the  applied  load  -  as  first  noted  by  Toya  (1974).  We  also  draw  attention  to  the  fact 
that  many  standard  tables  of  stress  intensity  factors  for  the  homogeneous  problem  ignore  this 
crucial  limitation. 

Mote  precisely,  we  present  a  slight  modification  of  the  England  (1966)  solution.  Of 
course,  this  solution  involves  overlapping  of  the  crack  faces  near  die  dps.  In  view  of  recent  work 
by  Rice  (1988),  He  and  Hutchinson  (1989)  and  Suo  and  Hutchinson  (1989)  we  endeavor  to  give  a 
proper  interpretation  of  the  solution  of  the  stated  problem.  But  we  remark  that  the  relationship 
between  the  present  work  and  a  properly  formulated  Comninou  (1977)  model,  allowing  for 
interfacial  contact  will  be  reported  at  a  later  date. 

The  paper  concludes  with  the  application  of  the  results  to  the  determination  of  the  reduction 
in  transverse  stiffness  of  a  dilute  unidirectional  ceramic  matrix  composite. 

Finally  we  remark  that  extensions  of  the  present  work  to  multiple  debond  cracks 
anisotropic  fibers  and  matrix  and  large  fiber  concentrations  are  possible  •  but  omitted  here. 

ANALYSIS 

We  consider  the  plane  strain  problem  of  an  infinite'  isotropic  matrix  with  moduli  Ki  ,p.i 
which  contains  a  circular  fiber,  of  radius  a,  with  moduli  Ki.Ri.  The  fiber  is  perfectly  bonded  to 
the  matrix  except  over  the  region  r»  a,  |0|  £  $.  The  system  is  loaded  by  uniaxial  tension  T  at 
infinity  at  angle  Q)  to  the  x-axis.  the  faces  of  the  crack  are  tractionfree. 

A  solution  of  the  problem  may  be  obtained  in  terms  of  the  Muskhelishvil  potentials 
$(z)  and  y(z).  Perhaps  the  simplest  form  is  obtained  if  we  work  with  the  potentials  $(z)  and  x(z) 
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where 


(1) 


X(z)-t«,(^)  +  V 


Let9t»Xi*$2*X2  denote  the  solutions  in  the  two  regions.  Then  continuity  conditions  over  the 
interface  imply  the  existence  of  two  functions  F(z),  G(z)  which  are  analytic  except  possibly  at  the 
origin,  infinity  and  over  the  cndc 

F(z)  -  ^j(z)  -  x^(z)  .  W  >  a  , 

F(z)  »^(z)-xt(z)  .  H<a  . 
and 

G(z)*^J-^1(z)  +  jJj-xl(2)  .  N>*  . 


G« -§♦»(*) +jfrx',«  .  H<»  • 

Moreover,  in  the  present  problem,  the  fact  that  the  creek  surface  is  traction  free  implies  that 

F(z)  -  A  +  B/z2  .  (6) 

With  the  help  of  the  preceding  definitions,  the  stated  problem  reduces  to  a  standard  Hilbert  problem 
for  G(z),  see  England  (1966).  who  also  gives  the  solution  in  closed  form. 

In  the  present  study,  quantities  of  particular  interest  are  the  asymptotic  neartip  stresses  and 
crack  opening  displacements.  It  is  not  difficult  to  show  that  near  the  tip  A  (see  Fig.  2.) 


(2) 

(3) 

W 

(5) 


F  (6,  to,  p)  In  sin  6 
■fSa 


_ r _ 

2a  sin  6 


C7) 


where 


a  _  ui  *  Kim 

Hi  +  K2Ml 


and 


(8) 

(9) 


°  •  (10) 

and  where  F  (6,  co,  p)  is  a  complicated,  but  known,  function.  Of  course  the  above  solution 
implies  overlapping  of  the  creek  faces,  and  thus  gives  an  estimate  of  that  pan  of  the  crack  length 
(rj  over  which  contact  take  place,  this  situation  has  been  carefully  and  extensively  discussed  by 
Rice  (1988)  in  the  case  of  an  interface  creek  between  two  half  spaces.  Guided  by  the  work  of  Rice 
(1988)  we  introduce  the  parameter 


2a0  (ID 

aniassenthat  small  scale  contact  at  a  particular  tip  will  occur  when  f  «  1.  In  such  circumstances 
Rice  (1988)  defines  a  classical  stress  intensity  factor  which  is  essentially  obtained  from  (7)  by 

deleting  the  tenn(r/2a  sin  For  similar  problems  He  and  Hutchinson  (1989)  and  Suo  and 
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Hutchinson  (1989)  have  argued  the  desirability  of  setting  0  *  0  in  equations  (7)  and  (8)  -  at  least  in 
the  case  of  two  half  spaces.  In  the  sequel  we  investigate  die  consequences  of  both  hypotheses. 

We  first  examine  the  range  of  parameters  for  which  the  interfacial  arc  crack  remains  open. 
Typical  results  are  shown  in  Fig.  1.  We  first  note  the  range  of  open  cracks  for  various  crack  sizes 
and  load  orientations  in  the  homogeneous  case.  For  a  given  crack  size  there  is  always  a  limitation 
on  the  orientation  of  the  applied  load  far  the  crack  to  remain  open.  In  the  inhomogeneous  problem, 
we  find  that  taking  f}  ■  0  implies  a  larger  range  of  open  cracks  than  we  obtain  for  small  scale 
contact  Further,  the  smaller  the  allowed  scale  of  contact  the  smaller  the  range  of  open  cracks. 

We  also  present  results  for  the  classical  stress  intensity  factors  for  a  debond  crack  at  the 
interface  of  an  SCS6  fiber  in  an  LAS  matrix  in  Figs.  2  and  3.  We  note  that  these  curves  apply  only 
when  the  cracks  are  open.  And  since  the  assumption  of  small  scale  contact  leads  to  a  smaller  range 
of  open  cracks  than  taking  0  =  0,  the  indicated  curves  in  Figs.  2  and  3  have  different  ranges  of 
validity.  However,  in  that  range  when  both  curves  are  valid,  there  is  little  between  the  S.LF.  in 
Figs.  2  and  3.  An  extensive  discussion  of  this  issue  will  be  given  elsewhere. 

Next  we  use  the  preceding  results  to  calculate  the  loss  of  transverse  stiffness  in  a  dilute 
unidirectional  SCS6/LAS  composite  due  to  interfacial  debonds.  For  simplicity  we  assume  that  the 
same  debond  crack  occurs  at  every  fiber  -  this  merely  has  the  effect  of  reducing  the  number  of 
geometrical  parameters.  In  addition  we  use  the  Budiansky  and  O'Connell  (1976)  definition  of 
crack  density: 

x 

where  N  is  the  number  of  cracks  per  unit  area  and  l  is  the  half  length  of  the  crack. 

It  must  be  emphasized  that  the  present  discussion  only  applies  to  open  cracks.  In  the  first 
instance  we  give  in  Fig.  4  the  results  for  a  homogeneous  material  -  so  that  the  overtapping  problem 
is  absent  Similar  results  are  presented  in  Figs.  5  and  6  for  the  composite  materiaL  Again  the 
difference  between  Figs.  5  and  6  merely  lies  in  the  range  of  validity  of  the  curves  for  which  the 
cracks  are  open. 

Finally  we  remark  that  an  extensive  discussion  of  the  issues  raised  herein  will  be  presented 
elsewhere. 
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FIBER  :  SCS  8  ;  E  3408  CP*,  v  3  0.23 
MATRIX  :  LAS  :  E  «83  CPfc  *  *  0 J 
fi  ■  0.0 


Fig.  3.  Classical  stress  intensity  factor  far  open  cracks  with  3-0. 


HOMOGENEOUS  MATERIAL 
CRACK  DENSITY,  c  =  0.03 


Fig.  4.  Loss  of  transverse  stiffness  u  a  function  of  orientation  of  applied 
load  for  open  cracks  in  homogeneous  material. 
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FIBER  :  SCS  0  :  E  *400  CP*,  v  *  0.23 
MATRIX  :  LAS  :  E  *03  CP*,  v  «  0.3 


r  *  o.oi .  Cf  *  o.i 


Load  orientation.  u 


Fig.  5. 


Los*  of  transverse  stiffness  as  *  funaion  of  orientation  of  applied 
load  for  open  ends  in  SCS6/LAS  and  small  scale  contact 


FIBER  :  SCS  9  ;  E  =400  GP*.  v  *  0.25 
MATRIX  :  US  ;  E  -03  CP*,  v  -  0.3 


d  a  IS  deg  rw< 
d  *  30  degree* 

#  s  45  degree* 

#  *  00  degree* 

•  •  it  ie.e  it.#  »•  »u  mi  stt  4#.f  til  sm  m  •  hi  k.o 
Load  orientation,  u 


Fig.  6.  Lon  of  transverse  stiffness  as  t  function  of  orientation  of  applied 
load  for  open  cracks  in  SCS&LAS  with  p  «  0. 
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